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It  has  long  been  suspected  that  cumulus  clouds  are  major  factors 
in  the  vertical  transport  of  energy,  mass,  and  momentum  in  the  atmos- 
phere, particularly  in  the  tropical  latitudes.  If  they  are,  those 
factors  must  be  taken  into  consideration  in  designing  simulation  models 
of  larger-scale  atmospheric  phenomena,  especially  those  depicting  the 
general  circulation  of  the  entire  world-wide  atmosphere.  Typically, 
planetary-scale  models  use  computation  grid  intervals  too  large  to  per- 
mit the  horizontal  resolution  of  cumulus  activity;  therefore,  the  cloud 
itself  cannot  be  directly  included  in  the  set  of  governing  equations. 

But  the  cumulative  effect  of  all  clouds  within  a grid  cell  must  be  taken 
into  account  if  the  larger-scale  model  is  to  properly  simulate  the  global 
atmospheric  dynamics  and  energetics.  For  that  reason  there  has  been  a 
considerable  effort  to  develop  parametric  schemes  to  account  for  cumulus 
ac  -vity , using  approaches  that  have  included  both  field  measurements 
and  theoretical  arguments. 

A useful  parametric  scheme  would  be  one  that  In  some  way  realisti- 
cally simulated  the  integrated  effect  of  the  whole  set  of  convective 
elements  on  the  large-scale  circulation,  temperature,  and  other  vari- 
ables. It  must  not  require  as  input  any  information  with  higher  resolu- 
tion in  space  or  time  than  the  standard  resolution  of  the  large-scale 
model,  nor  may  its  output  be  of  such  high  resolution.  If  the  parameter- 
ization were  simple  enough,  it  could  be  incorporated  into  the  large- 
scale  equations  and  solved  simultaneously  with  them.  More  probably,  it 
would  be  too  complex  for  that,  and  so  after  the  large-scale  computations 
were  made,  the  parameterization  would  take  these  preliminary  results  as 
input  data,  would  determine  the  changes  to  the  variables  due  to  convec- 
tion, and  would  feed  the  revised  conditions  back  to  the  large-scale  model. 

A reasonable  way  to  go  about  developing  such  a parameterization 
would  be  to  determine  first  in  some  detail  how  an  individual  convection 
element  affects  its  environment,  and  then  to  generalize  from  this  knowl- 
edge. The  present  work  represents  only  the  first  step  in  this  procedure. 
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In  this  work  the  transformation  from  one  kind  of  energy  (e.g., 
latent  enthalpy,  kinetic  energy)  to  another  and  the  rearrangement  in 
time  and  space  of  different  kinds  of  energy  as  a result  of  cumulus 
convection  are  analyzed  by  the  use  of  a numerical  model  of  cumulus 
dynamics.  These  kinds  of  data  provide  the  basic  information  necessary 
to  develop  parametric  schemes  for  use  in  the  larger-scale  models. 

This  report  is  one  of  a series  on  the  dynamics  of  climate,  a Rand 
project  sponsored  by  the  Defense  Advanced  Research  Projects  Agency. 

A large  part  of  the  work,  however,  was  supported  by  the  Navy  Environ- 
mental Prediction  Research  Facility  and  by  Corporate  funds  set  aside 
for  Rand-sponsored  restarch.  It  is  directly  related  to  Rand's  continu- 
ing interest  in  smaller-scale  atmospheric  phenomena.  For  other  related 
work,  see  the  following  Rand  publications: 

R-852-ARPA,  Numerical  Experiments  on  the  Relation  Between 

Microphysics  and  Dynamics  in  Cumulus  Convection, 

F.  W.  Murray  and  L.  R.  Koenig,  September  1971. 

R-877-ARPA,  A Documentation  of  the  Mintz-Arakawa  Two-Level 
Atmospheric  General  Circulation  Model,  W.  L. 

Gates,  F..  S.  Batten,  A.  B.  KahJe,  and  A.  B. 

Nelson,  December  1971. 

R-1005-ARPA,  The  January  Global  Climate  Simulated  by  the 
Two-Level  Mintz-Arakawa  Model:  A Comparison 

with  Observation,  W.  L.  Gates,  November  1972. 

R-1415-ARPA,  Some  Predicted  Climatic  Effects  of  a Simulated 

Sahara  Lake,  R.  R.  Rapp  and  M.  Warshaw,  March 
1974. 

R-1429-ARPA,  A Numerical  Experiment  on  the  Effects  of  Regional 
Atmospheric  Pollution  on  Global  Climate,  L.  R. 

Koenig,  June  1974. 
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SUMMARY 


Using  a numerical  model  previously  designed  to  simulate  atmospheric 
convection,  this  study  investigates  the  effects  of  cumulus  clouds  on 
the  temporal  and  spatial  characteristics  of  atmospheric  energy  content. 

The  forms  of  energy  investigated  are  latent  enthalpy,  thermal 
enthalpy  of  water  vapor,  thermal  enthalpy  of  dry  air,  thermal  enthalpy 
of  airborne  liquid  and  of  rain  that  has  fallen,  potential  energy  , nd 
kinetic  energy.  Work  done  by  the  pressure  force  is  excluded  because 
the  volume  of  the  domain  is  constant.  There  is  no  change  of  pressure 
at  the  upper  boundary,  and  only  very  small  changes  at  the  lower  and 
lateral  boundaries. 

Changes  in  energy  content  of  the  cloud  environment  and  the  flax 
of  energy  are  shown  as  functions  of  height  and  radial  distance  from 
the  cloud  axis.  These  changes  are  also  summed  vertically  and  hori- 
zontally over  the  domain  to  show  the  horizontal  and  verticil  rearrange- 
ment of  energy.  The  influence  of  the  cloud  is  largely  confined  to 
distances  within  a few  cloud  radii. 

The  relative  importance  of  different  energy  forms  is  a function 
of  position  with  respect  to  the  cloud.  Energy  related  to  the  presence 
of  water  vapor  (latent  plus  thermal  enthalpy  of  water  vapor)  has  been 
lound  to  account  for  the  majority  of  the  change  in  energy  in  the  vicin- 
ity of  the  cloud.  In  the  outer  region  changes  in  enthalpy  of  water 
vapor  are  counterbalanced  by  changes  in  the  enthalpy  of  dry  air,  but 
potential  energy  change  dominates. 

The  total  energy  used  in  attempting  to  decrease  conditional  in- 
stability in  the  atmosphere  is  given  by  the  total  amount  of  energy  pro- 
duced by  the  condensate  formed — regardless  of  whether  or  not  it  falls 
as  rain.  Clouds  of  differing  sizes  will  cause  differing  amounts  of 
energy  rearrangement.  The  simulated  "clouds"  used  in  this  work  re- 
semble in  size  the  cells  in  some  of  the  larger  natural  clouds  (but  not 
tK  largest  of  them)  observed  in  the  vicinity  of  the  station  at  which 
the  sounding  was  taken.  Thus  their  convective  activity  should  repre- 
sent reasonably  well  that  of  the  real  cloud  cells.  The  amount  by  which 
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the  action  of  the  computed  cell  reduces  the  potential  instability  of 
the  atmosphere  in  its  neighborhoou  is  0.5  percent  of  that  required  to 
establish  neutral  stability.  This  datum,  together  with  reasonable 
assumptions  about  the  number  and  frequency  of  occurrence  of  such  cells, 
leads  to  an  estimated  relaxation  time  of  convection  of  approximately 
ten  hours.  That  is,  ten  hours  is  required  for  the  departure  from 
neutral  stability  to  be  reduced  to  1/e  >f  its  initial  value. 
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(3-41) 

1 . INTRODUCTION 


Cumulus  convection  is  known  to  play  an  important  role  in  the  trans- 
formation and  rearrangement  of  energy  in  the  atmosphere.  As  used  here, 
"transformation"  means  the  change  of  one  type  of  energy  to  another 
(e.g.,  ootential  energy  to  kinetic  energy),  whereas  "rearrangement"  in- 
cludes both  changes  of  that  type  and  transport  from  one  location  to 
another.  The  conceptually  simple  transformation  of  latent  heat  of  con- 
densation to  sensible  heat  upon  the  condensation  of  v ater  vapor  grossly 
oversimplifies  the  true  cloud-environment  relationship.  In  actuality, 
the  release  of  latent  heat  causes  the  generation  of  kinetic  energy, 
which  In  turn  results  In  the  vertical  and  horizontal  advection  of  mass 
and  consequently  of  enthalpy  and  potential  energy.  There  are  th*n 
further  transformations  among  the  various  types  of  energy  as  a function 
of  space  and  time.  Presumably,  these  transformations  and  rearrangements 
.act  to  reduce  the  potential  instability  of  the  atmosphere.  One  of  our 
problems  is  to  find  out  how  this  is  accomplished. 

Although  the  amount  of  energy  involved  in  the  life  cycle  of  a 
single  convective  cell  is  negligible  by  comparison  with  the  energy  of 
the  larger  organized  circulations,  the  cumulative  effect  of  many  such 
cells  has  a profound  influence  on  these  larger  circulations.  Models 
of  the  general  circulation  or  even  mesoscale  models  cannot  resolve  con- 
vective cells,  and  so  they  must  parameterize  the  convective  activity  as 
a suogrid-scale  effect.  Such  a parameterization  must  realistically 
treat  the  cumulative  effect  of  an  ensemble  of  convective  cells  on  atmos- 
pheric energetics;  consequently,  it  should  be  based  on  the  knowledge  of 
the  energetics  of  a single  cell. 

Recognition  of  the  complexity  and  importance  of  cloud  energetics 
has  led  to  attempts  to  measure  in  the  field  those  factors  that  describe 
the  relationship  between  a cloud  and  its  environment  (e.g.,  Braham, 
1952),  and  to  a lesser  degree  to  use  theoretical  models  (e.g., 
Schlesinger,  1972)  to  analyze  the  cloud-environment  interactions  in 
terms  of  energy.  The  present  study  falls  in  the  latter  cate?  y. 
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This  study  makes  use  of  a numerical  model  designed  to  simulate 
the  life  cycle  of  an  individual  cumulus  cell  (Murray  and  Koenig,  1972). 
The  cloud  model  was  not  developed  for  the  study  of  energetics;  however, 


it  produces  fields  of  motion,  temperature,  pressure,  humidity,  and 
liquid  water  content  at  frequent  time  intervals,  which  makes  it  possible 
to  compute  the  content  of  energy  as  a function  of  time  and  space.  Hence, 
the  output  of  the  existing  cloud  model  was  used  as  the  input  for  the 
present  study  of  energetics. 

Briefly,  the  cloud  model  effects  the  numerical  solution  of  a set 
of  equations  including  the  equations  of  motion,  the  equations  of  con- 
servation of  mass  of  air  and  of  water  substance,  the  first  law  of 
thermodynamics,  the  equation  of  state,  and  several  equations  describ- 
ing microphysical  processes  related  to  condensation,  evaporation,  drop 
coalescence,  and  fallout  of  rain.  The  Boussinesq  approximation  (Spiegel 
and  Veronis,  1960)  is  used  to  simplify  the  hydrodynamical  and  thermo- 
dynamical equations,  and  the  parameterization  of  Kessler  (1969)  is  used 
for  the  microphysics.  The  solution  is  carried  out  in  a vertical  plane 
through  an  axis  about  which  horizontal  symmetry  is  assumed.  Initially, 
the  atmosphere  is  without  liquid  water,  is  horizontally  homogeneous,  and 
is  at  rest.  Convection  is  started  by  the  introduction  of  a small  per- 
turbation of  buoyancy  or  of  momentum.  Although  the  finite-difference 
scheme  was  not  designed  to  conserve  total  energy,  considerable  care  was 
taken  to  conserve  the  total  mass  of  water  substance,  which  consists  of 
vapor,  liquid  drops,  and  rain  that  has  fallen  to  the  ground.  Since 
water  substance  is  involved  in  much  of  the  cloud  energetics,  this  fea- 
ture contributes  strongly  toward  conservation  of  total  energy. 

"Tie  programs  designed  for  the  present  study  read  the  values  of  the 
meteorological  variables  computed  by  the  cloud  program  and  use  them  to 
compute  a number  of  other  quantities.  Among  these  are  specific  energy 
of  each  type  at  each  grid  point  and  specified  time,  changes  in  these 
values  from  time  zero  up  to  the  specified  time,  horizontal  and  vertical 
summations,  and  vertical  fluxes. 

Since  one  of  the  objectives  is  to  relate  cumulus-scale  convection 
to  circulations  of  mesoscale  through  planetary  scale,  some  attention 
must  be  paid  to  problems  of  resolution.  Atmospheric  circulations  of 
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horizontal  scale  larger  than  about  100  km  are  anisotropic  in  that 
their  motions  are  quasi-hor izontal . As  has  been  mentioned,  cumulus 
convection  is  completely  nonresolvable  in  the  horizontal  on  these 
larger  scales.  Their  vertical  and  time  scales,  however,  are  such  that 
cumulus  convection  can  be  crudely  resolved.  That  is,  vertical  trans- 
ports on  the  cumulus  scale  can  work  to  relieve  instabilities  in  the 
larger  scales,  and  the  length  of  time  required  for  a cumulus  cell  to 
go  through  its  life  cycle  is  comparable  to  the  appropriate  time  step 
for  the  larger  scales.  Hence,  this  study  emphasizes  average  vertical 
fluxes  and  changes  that  occur  from  the  time  of  inception  of  the  con- 
vective cell  until  the  conclusion  of  convection. 

The  term  "cloud"  is  frequeitly  used  in  this  work  to  designate 
the  airborne  condensate  simulated  by  numerical  computation.  Given 
that  the  model  simulates  the  initiation,  development,  and  decay  of  a 
convective  element  through  one  cycle  of  vertical  drafts,  the  term 
"cloud  tower"  may  more  correctly  describe  the  simulated  entity,  since 
in  nature  a cumulus  cloud  usually  consists  of  more  than  one  cell  or 
tower.  The  characteristics  of  a "cell"  simulated  herein  are  rather 
similar  to  those  of  a representative  thunderstorm  cell  as  described 
by  Brahar  (1952). 
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II.  CUMULUS  ENERGETICS 


ENERGY  CATEGORIES 


The  energy  content  of  a parcel  of  cloudy  air  may  be  subdivided 
into  various  categories: 


1. 

2. 


Kinetic  energy 
Static  energy 

a.  Potential  energy 

b.  Enthalpy 

(1)  Thermal  enthalpy 

(2)  Latent  enthalpy 


The  terms  "thermal  enthalpy"  and  "latent  enthalpy"  are  used  by 
van  Mieghem  (1973).  Meteorologists  ordinarily  use  the  terms  "sensible 
heat"  and  "latent  heat"  for  these  properties.  Thermal  enthalpy  is  re- 
lated to  the  heat  capacity  of  the  substance  involved,  and  latent  en- 
thalpy is  the  energy  related  to  a change  of  phase. 

In  thir  report  the  adjective  "specific"  applied  to  any  category 
of  energy  associated  with  a substance  will  indicate  the  energy  per 


unit  mass  of  that  substance.  Let  m,  be  the  mass  of  substance  d>  con- 


tained  within  a parcel  nf  cloudy  air,  where  "<j>"  may  be  replaced  by  "a" 
for  dry  air,  "v"  for  water  vapor,  or  "S."  for  condensed  water.  (See 
list  of  symbols  above.)  The  mixing  ratio  of  substance  <(>  is  defined  as 


m 


(2-1) 


If  h^  is  the  specific  energy  related  to  substance  <j),  then 


h* 


(2-2) 


represents  the  energy  of  that  substance  per  unit  mass  of  dry  air;  it 
will  be  designated  the  "dry-air  specific  energy."  Let  the  various 
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categories  of  energy  be  designated  by  subscripts,  as  follows: 

h gross  or  total  energy 

bj,  kinetic  energy 

hg  static  energy 
bp  potential  energy 
bp  enthalpy 
h^,  thermal  enthalpy 
hp  latent  enthalpy 

These  subscripts  maj  also  be  combined  with  the  subscripts  a,  v, 
and  l to  denote  the  substance  involved.  If  neither  a,  v,  nor  l is 
used,  the  symbol  refers  to  cloudy  air. 

Energy  expressed  in  the  form  of  (2-2),  being  referred  to  the  same 
unit  mass  for  all  categories,  is  additive.  Hence 


h*G  ' "*K  + "*S 

(2-3) 

h\.  - h*p  + h*E 

(2-4) 

h*  - h*  + h* 

E T L 

(2-5) 

h*E  = h*Ea  + h*Ev  + h*El 

(2-6) 

h*T  = h*Ta  + h*Tv  + h*T£ 

(2-7) 

Since  all  three  components  of  cloudy  air  share  the  same  height, 
and  except  for  fallout  of  rain,  the  same  velocity,  the  potential  and 
kinetic  energies  are  not  broken  down  by  substance.  The  enthalpy  of 
each  substance,  however,  is  computed  separately. 
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ENERGY  EQUATION  FOR_A 1_PARCEL ,_OF ^UWDY_AJR  IN  MOTION 

The  cloud  model  on  which  the  present  study  is  based  is  formulated 
in  terms  of  the  Boussinesq  approximation  as  set  forth  by  Spiegel  and 
Veronis  (1960),  but  with  modifications  to  allow  for  the  presence  of 
water  vapor  and  condensed  water  and  to  make  it  more  applicable  to 
moderately  deep  convection.  Spiegel  and  Veronis  divide  all  variables 
of  state  into  three  parts:  a mean  (constant)  part,  a part  dependent 

only  on  height  at  initial  time  in  the  absence  of  motion,  and  a part 
related  to  motion.  Thus,  for  density, 

0 = pm  + Pz(z)  + °'^x’  y»  z’  (2-8) 

The  equation  of  motion  then  would  be 

dv  , 

_=_l_Vp.  - gLj(  + vV2v  (2_9) 

m m 

where  v is  the  velocity  (having  horizontal  component  u and  vertical 
component  w),  and  the  other  symbols  are  as  defined  in  the  list  above. 

In  the  less  restrictive  anelastic  equations  of  Ogura  and  Phillips 
(1962),  the  reference  atmosphere  is  isentropic,  and  the  variables  of 
state  are  divided  into  only  two  parts: 

P = p(z)  • p' (x,  y,  z,  t)  (2-10) 

In  the  present  context,  p represents  the  density  of  cloudy  air;  hence 

p = p + p + p = p (q  + q + q ) (2-11) 

a v l a a Mv  MJr  v ' 

where  is  construed  as  a "bulk  density."  The  symbol  q has  been  re- 
tained  here  and  elsewhere  for  symmetry  and  clarity,  even  though  q El. 

3 

The  equation  of  motion  can  then  be  written  as 

d~  1 o’  2 

— = _ _ Vp'  - g k + vV^v 

P P 


(2-12) 
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\ 

The  advantage  of  (2-12)  over  (2-9)  is  its  better  capability  for 

deep  convection.  Also,  in  the  study  of  energetics  the  use  of  p gives 

better  results  than  does  the  use  of  p . The  cloud  model  used  herein 

m 

has  elements  of  both  systems  in  it,  but  (2-12)  will  be  used  to  develop 
the  energy  equation.  If  scalar  multiplication  by  v is  performed  on 
(2-12)  and  it  is  noted  that  v • k = w = dz/dt,  there  results 


d 

dt 


+ 


£± 

P 


V 


+ vv  • V v 


(2-13) 


Expansion  of  p in  the  manner  of  (2-10)  permits  the  first  term  on  the 
right-hand  side  of  (2-13)  to  be  written 


- ^ v • Vp ' = ~ “ (v  * Vp  - v • Vp) 


= Z («  £ - X * Vp) 


(2-14) 


Although  hydrostatic  balance  does  not  generally  hold  in  a convec- 
tive situation,  the  reference  atmosphere  is  at  rest  and  therefore  hydro- 
static. Hence  (2-14)  becomes: 

_ I y . Vp'  » - i (wgp  + V • Vp)  * - (gz)  - ^y  • Vp  (2-15) 

P P P 

Furthermore , 


1 „ 1 3p'  1 d£ 

- z X • ’p  - z at-  - z 

P P P 


(2-16) 


Combining  (2-13),  (2-15),  and  (2-16), 


d_ 

dt 


+ 


Pi 

P 


1 d£  1 V 

* - dt  - 3t 

P P 


+ vv  • 
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or 


d j + P + P1  d 

dt  \ 2 / - dt 

\ ' p 


(gz)  + 


I 


dt  - 3t 

P 


vv 


V2v 


(2-17) 


The  terms  on  the  left-hand  side  of  (2-17)  will  be  recognized  re- 
spectively as  the  rate  of  change  of  kinetic  energy,  potential  energy, 
and  enthalpy  (provided  there  are  no  external  sources  or  sinks  of  heat, 
such  as  radiation).  All  refer  to  unit  mass  of  cloudy  air.  The  terms 
on  the  right-hand  side  represent  the  work  performed  on  the  parcel  by 
the  pressure  force  and  the  eddy  viscosity  force.  Both  are  small. 

Since  the  third  term  represents  the  rate  of  change  of  total  en- 
thalpy of  a unit  mass  of  cloudy  air,  it  can  be  written  dh^/dt.  It 

remains  to  evaluate  h„. 

E 

In  accordance  with  (2-6)  together  with  (2-1)  and  (2-2),  the  en- 
thalpy of  a parcel  of  cloudy  air  is 


(m  + m + m )h 
a v HE 


m h_  + 
a Ea 


mvhEv  + m*hE£ 


(2-18) 


The  mass  of  the  air  remains  constant,  and  except  for  the  fallout 
of  rain,  which  for  the  time  being  will  be  neglected,  the  total  mass  of 
water  substance  remains  constant.  That  is, 

dm  = -dm.  (2-19) 

v l 


Hence,  differentiation  of  (2-18)  leads  to 


(m  + m + m„)  dh_  = m dh„  + m dh„  + m dh„„ 
a v l E a Ea  v Ev  l Ef. 


+ (hEv  - hEi)  dmv 


(2-20) 


Expanding  (2-20)  according  to  the  cnain  rule  of  dif ferentiati  ->n. 
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(m  + m + ra  ) dh  = m (— -|ja) 
a v £ E a \ 3T  / 

dT  + "a(  3p 

+ m 

V 

f9hEv' 

l 9p  t 

) dp 

+ Bl 

' ahE«  ) 

i Sp  ) 

1 dp 

+ (h  - h ) dm 

Ev  E£  v 

(2-21) 

The  changes  in  specific  enthalpy  at  constant  pressure  are  by 
definition  the  specific  heat  capacities: 


The  changes  in  specific  enthalpy  at  constant  temperature  are  zero, 
and  the  difference  in  specific  enthalpy  between  the  two  phases  is  the 
latent  heat  of  condensation: 


Hence  (2-21)  becomes 


(m  + m + m„)  dh„ 
a v Jc  E 


me  dT  + m c 
a pa  v pv 


dT 


+ me  dT  + L dm 

x.  Jc  v 


(2-23) 


12-24) 


or,  after  division  by  m 

a* 

dh*F  = (qa  + qv  + q£>  dhE  = (qaCpa  + Vpv  f Vi>  dT  + L dqv  (2‘25) 

Following  the  method  of  Iribarne  and  Godson  (1973),  we  can  use 
(2-23)  to  transform  (2-18)  into 
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(ma  + \ + m£)hE 


m h + ra  L + (m  + m.)h„„ 
a La  v v £ El 


(2-26) 


Differentiation  of  (2-26)  In  the  fashion  of  (2-21),  followed  by 
division  by  m^,  yields 

dh*E  = (qa  + qv  + q^>  dhE 

" (qaCpa  + (qv  + qf)c£)  dT  + L d%  + qv  dL  (2"27) 

The  two  expressions  for  the  dry-air  specific  enthalpy  of  a parcel 
of  cloady  air,  (2-25)  and  (2-27),  are  equivalent.  Ordinarily,  the 
variation  of  latent  heat  with  temperature  can  be  neglected  in  meteor- 
ology,  but  if  that  is  done  in  integrating  (2-27),  the  resulting  ex- 
pression for  the  dry-air  specific  enthalpy  of  water  vapor  is  q c„T, 

^v  £ 

whereas  the  correct  expression  from  the  integration  of  (2-25)  is 
qvCpvT*  The  former  ls  to°  large  by  a factor  of  cx/cpv  = 4187/1846 
- 2.27.  If  (2-25)  and  (2-27)  are  equated,  Kirchoff's  equation  results: 


c 


£ 


(2-28) 


In  order  to  make  (2-25)  compatible  with  (2-17),  ii.<?  enthalpy  must 
be  expressed  in  terms  of  unit  mass  of  cloudy  air  rather  than  of  dry  air. 
This  can  be  accomplished  by  multiplying  (2-25)  by  p /p  or  alternatively 
by  dividing  (2-24)  by  m + m + m..  In  accordance  with  (2-11)  and  the 

<1  V X, 

conservation  of  mass  of  air  and  of  total  water  substance. 


qa  + qv  + q £ 


Hence  from  (2-25) 


(2-29) 


P P 

dh  = — (qc  +qc  + q c ) dT  + ■—  L dq 
E — a pa  v pv  £ £ — Mv 


(2-30) 
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and  since 


(2-31) 


it  follows  that  (2-17)  becomes 


d 

dt 


P d , . a 

- k <8*>  + — 


(q  c 
a pa 


+ 


q c 
v pv 


+ Vt> 


dT 

dt 


1 3El 

- 3 1 
P 


+ vv 


(2-32) 


The  terms  on  the  left-hand  side  of  (2-32)  represent  the  rate  of 
change  of  foui  kinds  of  specific  energy:  kinetic  energy,  potential 

energy,  thermal  enthalpy,  and  latent  enthalpy.  Therefore,  using  (2-2) 
and  (2-29) 


h*K  “ (qa  + % + <Ig)(u2  + w2)/2  (2-33) 

h*P  = (qa  + qv  + qH)  - gz  (2-34) 

P 


h*  = 
T 

(q  c + q c + q.c.)T 

a pa  nv  pv  M£  £ 

(2-35) 

h*Ta 

q c T 
a pa 

(2-36) 

h*Tv  " 

q c T 
v pv 

(2-37) 

h*T£ 

ViT 

(2-38) 

h*L  = 

%L 

(2-39) 

The  factor  p/p  in  (2-34)  is,  in  accordance  with  (2-29),  equal  to 
Pa/pa’  which  is  very  nearly  unity,  so  it  might  well  be  neglected.  As 
will  be  seen  later,  however,  the  nature  of  the  Boussinesq  approxima- 
tion dictates  that  it  should  be  retained. 
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All  forms  of  energy  shown  In  (2-17)  with  the  exception  of  those 
related  to  local  pressure  change  and  eddy  viscosity  are  defined  by 
(2-33)  through  (2-39).  In  a Boussinesq  system  the  local  pressure 
change  is  small  and  is  usually  ignored.  The  cloud  model  computations 
confirm  that  it  is  extremely  small,  and  since  the  highly  implicit 
nature  of  the  problem  makes  that  computation  less  accurate  than  those 
of  the  other  variables  of  state,  the  local  pressure  change  term  will 
not  be  further  discussed.  The  eddy  diffusion  term  presumably  repre- 
sents a conversion  of  kinetic  energy  to  heat,  but  in  a finite-difference 
solution  it  could  represent  a sink  of  energy  through  the  cascade  effect 
to  larger  wave  numbers.  In  any  case,  the  kinetic  energy  itself  is  very 
small,  so  the  eddy  diffusion  is  quite  negligible. 


ENERGY  DENSITY 

For  purposes  of  summation,  it  is  more  convenient  to  refer  energy 
to  unit  volume  rather  than  to  unit  mass.  This  can  be  achieved  by 
multiplying  (2-32)  by  p,  resulting  in 


P (q  c 

a a _>r 


+ q c + q c„) 
v pv  M l V 


d_T 

dt 


dq  » 

+ pL^v  = |2  + 


dt 


at 


vpv 


V2v 


(2-40) 


Then,  corresponding  to  the  expressions  for  energy  per  unit  mass  of  dry 
air  found  in  (2-33)  through  (2—39)  , one  may  write  the  following  ex- 
pressions for  energy  per  unit  volume,  or  energy  density,  with  the  aid 
of  (2-11)  and  (2-29): 

\ = Pa(qa  + qv  + q£)(u2  + \2)/ 2 (2-Ai) 


np  = Pgz  = Pa(qa  + qv  + q£)gz 


(2-42) 


hT  15  p (q  c + q c + q.c.)T 
T a a pa  v pv  Hi 


(2-43) 
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nTo  “ 0 q c T 
Ta  a a pa 


nT„  = p q c T 
Tv  a v pv 


nTt  ’ °a'1lcilT 


n.  » D Lq 
L a v 


(2-A4) 

(2-45) 

(2-46) 


(2-47) 
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IJLL_  COMPUTATIONAL  METHO  _S 


VOLUME  ELEMENTS 

The  Rand  cloud  model  (Murray  and  Koenig,  1972)  as  used  in  the 
present  experiment  is  formulated  in  cylindrical  coordinates  and  assumes 
axial  symmetry.  A grid  of  uniform  mesh  length  D is  constructed,  with 
the  horizontal  index  i running  from  1 at  the  central  axis  to  1 at  the 
lateral  wa.l,  and  the  vertical  index  k running  from  1 at  the  top  to  K 
at  the  bottom.  In  the  present  experiment,  D = 200  m,  I = 31,  and 

K = 45;  hence,  the  computational  domain  is  a cylinder  of  radius  6 km 
and  height  9 km. 

It  is  assumed  hat  values  for  3 grid  point  are  representative  of 
a volume  element  about  the  grid  point  delimited  by  vertical  cylindrical 
and  Horizontal  plane  surfaces  constructed  midway  between  the  given  grid 
point  and  its  closest  neighbors.  Thus,  the  volume  element  represented 
by  grid  point  (i,k)  has  the  horizontal  area 

= ttD^/4 

2 2 
Ai  * *(ri  + §' > * ~ y) 

* 2TrriD;  1 < i < I 

Ai  ■ "ri2  - ’(ri  - f> 

- »D(rj  + £) 

and  depth 

\ - D;  1 < k < k] 


D/2 


(3-2) 


r 
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where  r^,  the  radial  distance  irom  the  central  axis,  is  given  by 

= (i  - DD  (3-3) 

The  volume  of  the  element  represented  by  grid  point  (i,k)  is  then 


The  special  cases  for  i = I , k = 1 , and  k = K are  the  result  of  assum- 
ing  that  the  computational  domain  ends  at  rigid  surfaces  constructed 
through  th*  grid  points  on  the  boundary  rather  than  at  a distance  of 
D/2  beyond  them.  The  choice  of  this  boundary  assumption  was  arbitrary. 

SUMMATIONS 

The  total  energy  of  any  given  type,  0,  residing  within  the  volume 
element  represented  by  the  g-id  point  (i,k)  is  then 


i.k 


e.  , i,k 
i,k 


(3-5) 


where  0 can  be  replaced  with  K for  kinetic  energy,  P for  potential 
energy,  etc.  These  values,  then,  can  be  summed  in  various  ways.  In 
particular,  three  types  of  summation  are  used  in  the  present  work. 

Let  <J>  be  any  quantity  defined  at  the  grid  points.  Then  the  sum- 
mation of  (J)  by  stubs  (i.e.,  horizontally)  will  be  denoted  by  angular 
brackets : 


I 

i=l 

The  summation  of  <£  by  tubes  (i.e.,  vertically)  will  be  denoted  by 
square  brackets: 

K 

■♦'i  -12  *i,k 

k=l 


(3-7) 
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The  summation  over  the  entire  domain  will  be  denoted  by  braces: 

I K 

M ‘Y.  Y b.k  °-8) 

i=l  k=l 

It  will  be  noted  that  summation  over  the  domain  is  equivalent  to 
successive  summation  by  slabs  and  by  tubes  in  either  order;  i.e., 

j*|  = = [<«)k]  (3-9) 

If  the  total  mass  of  dry  air  in  the  closed  computational  domain 
is  to  be  conserved,  then  it  is  necessary  that 

I PaV < = constant  (3-10) 


This  is  assured  in  the  present  model,  which  is  based  on  the  Boussinesq 
approximation,  by  substituting  p , a i notion  of  k only,  for  p . This 

H cl 

substitution  is  fully  consistent  with  the  use  of  p in  computing  all 

Si 

energy  densities  except  potential  energy  density;  see  equation--  (2-41) 
through  (2-47).  It  imposes  the  restriction  that  the  work  of  expansion 
is  zero,  which  is  consistent  with  the  fixed  volume  of  the  domain  and 
the  fixed  volume  represented  by  each  grid  point. 

ENTHALPY  OF  RAIN 

During  the  operation  of  the  cloud  model,  a part  of  the  condensed 
liquid  falls  out  as  rain.  This  represents  a certain  amount  of  energy 
(thermal  enthalpy)  that  is  removed  from  the  atmosphere. 

If  is  the  depth  of  rain  on  the  ground  in  millimeters  over  the 

area  A.,  then  R.A  is  the  mass  in  kilograms  of  that  amount  of  water, 

1 11  3-3 

since  the  density  of  water  is  10  kg  m . The  thermal  enthalpy  re- 
moved from  the  atmosphere  by  rain  during  any  given  time  step  depends 
on  the  temperature  of  the  rain  during  that  time  step;  hence  the  enthalpy 
removed  from  tube  i between  time  0 and  time  x is 
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T 


C«(4Ri,t,AiTl,K,t 


(3-11) 


where  ( A R , ) is  the  amc  int  of  rain  falling  during  the  time  step  end- 

ing  at  time  t.  The  appropriate  rain  temperature  is  assumed  to  be  the 
air  temperature  in  the  lowest  level  at  that  time.  (Throughout  this 
study,  it  is  assumed  that  all  the  constituents  of  a parcel  of  cloudy 
air  at  a given  position  and  time  have  the  same  temperature.) 

Since  the  enthalpy  of  rain  is  not  a function  of  k,  summation  by 
tubes  is  a redundant  operation 


(3-12) 


Summation  by  slabs  is  equivalent  to  total  summation: 


<V 


(3-13) 


ENERGY  CHANGE 

Although  the  distribution  of  energy  by  type  and  location  at  a 
given  time  is  of  interest,  the  changes  in  the  course  of  the  cloud  sim- 
ulation are  of  greater  interst.  Let  the  operator  A represent  change 
from  time  zero  up  to  some  given  time  T,  or 

A$  - $t  - *0  (3-14) 

This  operator  can  be  applied  to  energy  densities,  as  given  in 
(2-41)  through  (2-47),  the  corresponding  energy  totals,  as  indicated 
in  (3-5),  or  energy  summations,  as  indicated  in  (3-6)  through  (3-8). 
Since  the  operator  A is  independent  of  the  operators  ( ),  [ ] , and 
{ },  the  order  of  application  is  immaterial. 

In  some  instances,  however,  straightforward  application  of  (3-14) 
leads  to  the  computation  of  small  differences  of  large  quantities,  with 
consequent  loss  of  precision.  Wherever  possible,  this  problem  was 
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minimized  by  factoring  out  those  parts  of  <J>  that  are  not  time-dependent. 
This  was  facilitated  by  the  decomposition  of  the  variables  of  state  in 
accordance  with  (2-10).  For  instance,  T tends  to  be  two  orders  of  mag- 
nitude larger  than  T',  so  it  is  more  accurate  to  compute  AT'  than  AT, 
even  though  they  are  ostensibly  identical. 

The  reference  values  of  u,  w,  and  q^  are  all  zero,  and  their  per- 
turbed values  are  small.  In  earlier  versions  of  the  model,  q was 

v 

treated  as  though  its  reference  value  were  zero  even  though  q^  varies 
in  a narrow  range  about  a value  some  two  orders  of  magnitude  larger 
than  q^.  In  a later  version,  q^  was  decomposed  in  the  same  manner  as 
the  other  variables,  but  this  appeared  to  have  little  effect  on  the 
energy  changes. 

In  practice,  application  of  (3—14)  to  (2-41)  through  (2-47)  was 
carried  out  as  follows: 

A\  * \ p£.  + % + + w2))  (3-15) 


(3-16) 
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After  (3-15)  through  (3-20)  have  been  multiplied  by  V to  derive 
. _ 1 }k 

, • ®tc.,  they  can  be  suoned  by  slabs  and  tubes. 

1 f K 

In  the  present  study,  rainfall  was  always  zero  at  time  zero,  so 


(3-21) 


FLUX 

If  nQ  is  the  energy  density  of  some  particular  category  0,  then 
the  corresponding  vertical  flux  of  energy  is 

Fe  - wne  (3-24) 

Any  variable  property  $ can  be  divided  into  a time-averaged  part 
and  a departure: 


(3-25) 

where,  if  there  have  been  n - 1 time  steps  between  the  initial  time 
and  the  final  time  T, 


$ 


(3-26) 


It  should  be  noted  that  the  decomposition  indicated  by  (3-25)  is  not 
the  same  as  that  indicated  by  (2-10) . The  straight  bar  (“)  implies 
a horizontal  space  average  at  initial  time,  whereas  the  tilde  O repre- 
sents a time  average  at  a given  location. 

'rom  (3-25)  and  (3-26)  it  is  readily  demonstrated  that 


♦ - ♦ 

♦"  - 0 
a$  ■ a$; 


(3-27) 


a Independent  of  t 
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The  flux  given  by  (3-24)  for  a specified  location  (i,k)  and  time  t 
can  be  decomposed  :'n  accordance  with  (3-25) 


or 


V.,'  (“  + “"4,k("a+  Vt) 

'i,k 

)iiM-(^)ijk+(VJik+(wVe)ijk+(wV)i>tt 


(3-28) 


If  (3-28)  i9  time-averaged  in  accordance  with  (3-26),  and  (3-27) 
is  used,  there  results 


) = (w0  + (w"n9”) 
}i,k  ' 0/i,k  ' 9 'i,\ 


(3-29) 


Since  n0  is  an  energy  density,  the  dimensions  of  (3-29)  are  J m-2  s"1. 

The  terms  cn  the  right-hand  side  of  (3-29)  can  be  identified  as  the 
flux  of  the  mean  energy  of  type  0 with  the  mean  vertical  wind  (or,  more 
simply,  the  mean  vertical  flux)  and  the  eddy  vertical  flux.  The  latter 
expresses  the  correlation  between  the  vertical  wind  and  energy  of  type  0. 

If  each  term  of  (3-29)  is  multiplied  by  A^,  the  product  is  the 
vertical  flux  through  the  entire  horizontal  surface  represented  by  the 
grid  point  (i,k).  Horizontal  summation  then  gives  the  total  flux 
through  the  level  designated  by  k: 


<AF0>  = <w  AnQ>  + <w"  An  ">  (3-30) 

k k k 

It  is  convenient  to  define  two  types  of  horizontal  mean.  The 

ordinary  one  fot  any  function  $ of  i and  k is 

1 , k 


(3-31) 
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sinot*  the  number  of  terms  in  die  slab  summation  is  I.  The  area-weighted 
horizontal  mean  is 


(3-321 


It  follows  that 


(3-33) 


r » 

(Note  that  4^  and  may  also  be  functions  of  time,  but  <j>  represents 
the  reference  state  at  initial  time.)  By  analogy  to  (3-25),  any  func- 
tion of  i can  be  broken  into  a slab-averaged  part  and  a departure: 


<t> 


i,k 


(3-34) 


If  we  make  the  definition 


h ‘ Ane 


(3-35) 


and  note  that  A is  independent  of  t,  (3-30)  becomes,  through  the  use 
of  (3-31), 


I (AF  ) = I (w£  ) + I (w"F,  •') 


(3-36) 


By  means  of  (3-34) 
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With  the  use  of  identities  analogous  to  those  of  (3-27),  the  slab  aver- 
age of  (3-37)  becomes 


■ \ \ + <S‘C>[ 


(3-38) 


Substitution  of  (3-38)  into  (3-36),  followed  by  use  of  (3-33),  yields 


KAFJ  -w 
k 


I An.  + I A(w*  n *)  ♦ I A(w"  n ") 

k 0 k 0 k 


(3-39) 


If  (3-39)  ^s  divided  by  (A)  and  (3-32)  is  used,  there  results 


~ “ ~ “ ~ ^ — < “ 

uk  \ + <”  "o  > + V> 

k k k 


(3-40) 


In  (3-40)  Fg  represents  the  vertical  flux  of  energy  of  type  0 
k 

through  the  horizontal  surface  designated  by  k,  averaged  both  hori- 

—2  —1 

zontally  and  in  time.  The  dimension  is  J m s , and  it  yields  a 
single  value  for  each  level  k.  The  first  term  on  the  right-hand  side 
is  the  flux  of  the  doubly  averaged  energy  by  the  doubly  averaged  verti- 
cal wind;  it  can  be  called  the  mean  flux  over  the  slab  as  distinguished 
from  the  corresponding  term  of  (3-29),  which  is  the  mean  flux  at  a grid 
point.  The  other  two  terms  of  (3-40)  represent  two  types  of  eddy  flux. 
The  first  of  these  two  is  associated  with  standing  waves  (or  circula- 
tions approximating  them).  Inspection  of  Figs.  13  and  14  reveals  that 
such  circulations  do  exist  in  this  model.  The  remaining  term  represents 
transient  eddies.  . 

The  summation  Z Fq  can  be  thought  of  as  representing  the  rise 
k=l  k 

or  fall  of  the  "center  of  energy"  (analogous  to  the  center  of  mass)  of 

energy  of  type  0.  Some  computations  of  this  quantity  are  shown  in 
Fig.  20. 


I 
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The  entire  foregoing  development  for  energy  flux  is  also  appli- 
cable to  mass  flux.  For  example,  to  study  the  flux  of  water  vapor,  it 

is  merely  necessary  to  replace  r|  in  the  above  equatic  s with  i q . 

a v 

The  accumulation  or  depletion  of  energy  in  a g^ven  layer  is  found 
by  vertical  differencing  of  (3-30).  This  is  a measure  of  the  divergence 
of  energy,  and  it  may  be  represented  by 

■WV  ' <Ah>  - <Ah>  o-«n 

k k+1 

Since  the  indexing  system  is  oriented  so  that  k increases  as  z 
decreases,  (3-41)  will  have,  a negative  value  in  the  case  of  convergence 
and  a positive  value  in  the  case  of  divergence.  That  is,  i:  (3-41)  is 
negative,  it  indicates  an  increase  in  energy  of  type  6 in  the  layer 
between  k and  k+1.  Of  course,  none  of  the  flux  calculations  tell 
anything  about  any  transformations  of  energy  that  might  be  occurring. 
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IV.  ENERGY  DISTRIBUTION  AND  TRANSFORMATION 

Materials  used  in  the  present  analysis  comes  from  two  cloud- 
simulating  computations  reflecting  different  stages  in  development  of 
the  study.  Table  1 summarizes  the  initial  conditions,  some  properties 
of  the  clouds  produced,  and  certain  aspects  of  the  computation.  Cloud 
B is  computed  from  the  San  Juan  sounding  of  20  August  1963  (see  Fig.  1) 
and  is  similar  to  Cloud  3 of  Murray  and  Koenig  (1972).  Cloud  A differs 

Table  1 


SOME  CLOUD  PROPERTIES 


Property 

Cloud  A 

Cloud  B 

Atmospheric  sounding 

Smoothed,  San  Juan, 
2300  GMT,  20  Aug.  1963 

Unsmoothed,  San  Juan, 
2300  GMT,  20  Aug.  1963 

Perturbation 

Momentum 

Humidity  (net  change 
equal  to  zero) 

Base 

1.2  km 

1.2  km 

Height 

Maximum 

Time  of  maximum 

5.8  km 
25  min 

6.0  km 
35  min 

Radius 

Maximum 

Time  of  maximum 

0.6  km 
20  min 

0.6  km 
30  min 

Temperature  departure 
Maximum 

Time  of  maximum 

2 . 1 °C 
15  min 

2 . 3°C 
25  min 

Vertical  draft 
Maximum 

Time  of  maximum 
Height  of  maximum 

18  m sec  * 
20  min 
A . 2 km 

i ^ -1 

15  m sec 
30  min 
4.8  km 

Life  time 

50  min 

55  min 

Total  condensate 

3558  tons 

Rain  amount 

282  tons 

4693  tons 

Conservation  applied? 

Yes 

No 

Eddy  diffusion? 

Yes 

Yes 

NOTE:  Time  step  £ 12  sec 


Time  resolution  = 5 min 
Distance  resolution  = 0.2  km 
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from  Cloud  B In  three  main  characteristics:  (a)  the  sounding  was 

slightly  smoothed  with  respect  to  temperature  and  strongly  smoothed 
with  respect  to  dewpoint  (see  Fig.  2);  (b)  the  initial  perturbation 
was  of  momentum  rather  than  humidity;  and  (r)  a computational  scheme 
was  included  to  assure  the  conservation  of  mass  of  water  substance. 
Because  of  this  third  character i st ic , Cloud  A is  considered  to  be  the 
more  useful  of  the  two.  Cloud  B is  used  mainly  to  study  the  vertical 
flux  of  various  categories  of  energy,  a feature  that  is  not  highly 
sensitive  to  either  conservation  of  mass  of  water  substance  or  ncn- 
smoothed  soundings  but  that  is  sensitive  to  a perturbation  of  the 
vertical  wind. 


Fig  1-Sounding  for  San  Juan,  2300  GMT,  20  August  1963.  used  for  Cloud  B 
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Fig.  2 Smoothed  version  of  San  Juan  sounding,  used  for  Cloud  A 


The  initial  impulse  of  Cloud  A provided  a circulation  with  a max- 
imum updraft  of  3.5  m sec  1,  thereby  introducing  7.5  * 109  joules 
(0.0075  terajoules)  of  kinetic  energy  into  the  computational  domain. 

If  no  condensation  occurs,  an  impulse  such  as  this  rapidly  decays  and 
there  is  little  rearrangement  of  the  various  forms  of  energy.  The 
humidity  impulse  of  Cloud  B introduced  more  energy  in  the  form  of 
enthalpy  into  the  active  cloud  region,  but  this  was  counterbalanced 
by  an  equal  reduction  in  enthalpy  in  the  environment.  The  momer turn 
impulse  of  Cloud  A caused  a more  rapid  cloud  development  than  did  the 
humidity  impulse  of  Cloud  B. 

The  evolution  of  Cloud  A is  summarized  in  Table  2.  The  top  of 
the  table  depicts,  at  five-minute  intervals,  the  total  mass  of  condensed 
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water  in  the  domain,  both  airborne  (cloud  water  mass)  and  on  the  ground 
(rain  mass),  and  also  the  cumulative  mass  of  water  that  has  been  con- 
densed up  to  the  given  time.  This  is  always  larger  than  the  sum  of  the 
other  two  categories  because  some  of  the  condensate  later  evaporates. 

The  lower  part  of  Table  2 shows  the  change  in  energy  content  of  various 
types  from  initial  time  up  to  15,  30,  45,  and  60  minutes.  The  right- 
hand  column  gives  the  base  level  from  which  the  changes  occur.  Except 
for  the  thermal  enthalpy  of  liquid  and  the  kinetic  energy,  the  base 
levels  are  all  much  larger  than  the  changes.  The  degree  to  which  energy 
is  conserved  is  shown  by  the  totals  to  be  1 part  in  10^. 

Essentially  no  new  condensation  occurred  after  30  minutes;  i.e., 
the  cloud  was  virtually  dead  from  that  time  on.  The  latent  enthalpy, 
thermal  enthalpy,  and  potential  energy  changes  also  indicate  that  the 
cloud  is  in  decline  after  30  minutes.  The  kinetic  energy  is  nearly 
constant  during  the  later  period,  but  it  is  less  than  the  initial  per- 
turbation, and  examination  of  the  fields  of  motion  shows  that  the 
circulations  are  much  less  organized.  Rainout  and  evaporation  erode 
the  cloud  until  it  finally  disappears  after  55  minutes. 

The  rearrangements  of  energy  of  Cloud  A are  shown  in  Figs.  3 through 

7 in  the  form  of  fields  of  change  of  energy  density,  expressed  in  joules 

per  cubic  meter  per  hour.  As  shown  in  Table  ] , the  volume  actually 

occupied  by  the  cloud  extends  from  1.2  km  up  to  6 km  and  from  the  axis 

out  to  0.6  km.  These  limits  approximate  the  envelope  of  the  cloud 

boundary  over  the  lifetime  of  the  cloud.  Figure  8 shows  the  change  of 

energy  content  of  each  type  after  being  summed  over  horizontal  slabs 

in  accordance  with  (3-6)  , and  Fig.  9 shows  the  corresponding  changes 

after  summation  over  vertical  tubes  in  accordance  with  (3-7).  When 

12 

summed  by  slabs  or  tubes,  energy  is  expressed  in  terajoules  (10  joules) 
and  energy  change  in  terajoules  per  hour.  All  changes  shown  in  these 
figures  are  taken  from  time  zero  (before  the  inception  of  convection 
but  after  the  implementation  of  the  perturbation)  to  60  minutes  (after 
completion  of  the  life  cycle  of  the  cell) . 

LATENT  ENTHALPY 

Changes  in  the  latent  enthalpy  are  shown  in  Figs.  3,  8(a),  and  9(a) 
for  Cloud  A.  Increases  in  latent  enthalpy  occur  in  regions  that  become 
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Fig.  5 -Change  from  0 to  60  min  of  enthalpy  density  Fig.  6-  Change  from  0 to  60  min  of  potential  energy  density 

of  dry  air  (J  m ’ hr  ' ).  Cloud  A (j  m 3 hr-' ),  Cloud  A 


9 


Radial  distance  (km) 

NOTE  Regions  of  decrease  are  shaded 

Fig  7— Change  from  0 to  60  min  of  static  energy  density 
IJm’  hr  1 ),  Cloud  A 
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Radiol  distance  (km) 

Fig.  9 — (e  - h ) 

richer  in  water  vapor  at  the  expense  of  those  that  become  poorer,  for 
latent  enthalpy  is  directly  related  to  the  mass  of  water  vapor.  As 
revealed  in  Fig.  3,  the  entire  region  beneath  the  cloud,  and  a large 
volume  encircling  the  cloud  and  extending  vertically  from  the  cloud 
base  to  about  half  the  height  of  the  cloud,  has  decreasing  latent  en- 
thalpy, the  result  of  downward  advection  of  drier  air  in  the  region 
distant  from  the  cloud  and  of  upward  advection  of  moister  air  into  the 
cloud  itself.  The  dying  cloud  leaves  behind  a pocket  of  increased 
latent  enthalpy  that  is  larger  in  volume  than  the  cloud  itself.  Near 
the  cloud,  the  changes  in  latent  enthalpy  are  relatively  large  in  mag- 
nitude (maximum  values  near  4000  J m"3),  and  the  gradients  are  intense 
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The  overall  vertical  movement  of  latent  enthalpy  is  shown  in 
Fig.  8(a).  There  are  increases  in  latent  enthalpy  at  levels  below 
that  of  the  cloud  base  (~  1.2  km)  and  throughout  the  volume  in  which 
the  upper  half  of  the  cloud  resides.  These  are  offset  by  decreases 
through  the  volume  in  which  the  lower  half  of  the  cloud  resides.  (The 
slabs  over  which  enthalpy  is  summed  cover  the  entire  horizontal  extent 
of  the  domain,  not  just  the  part  occupied  by  the  cloud.)  'T'he  subcloud 
increase  is  in  the  region  several  cloud  radii  away  from  the  cloud  it- 
self and  apparently  is  caused  by  convergence  of  the  horizontal  outflow 
associated  with  the  moist  downdraft.  The  middle  level  decrease  occurs 
in  the  region  on  which  the  cloud  feeds,  and  the  upper  level  increase 
is  caused  by  the  evaporation  of  water  in  the  dying  cloud.  The  net 
change  is  negative  because  of  the  removal  of  rainwater  formed  by  con- 
densing water  vapor. 

The  overall  horizontal  movement  of  latent  enthalpy  is  shown  in 
Fig.  9(a).  These  data  show  a marked  increase  in  latent  enthalpy  in 
the  cylinder  occupied  by  the  cloud  coupled  with  decreases  in  the  air 
away  from  the  cloud.  (The  downturn  in  all  the  curves  of  Fig.  9 as  the 
central  axis  is  approached  is  a result  of  the  dependence  on  the  square 
of  radial  distance;  see  Eq.  (3-1).)  About  2 cloud  radii  from  the 
central  axis  (i.e.,  about  1.3  km)  there  is  a minimum  point  on  the  curve, 
and  there  is  a relative  maximum  at  3 cloud  radii  (2.3  km).  Reference 
to  Fig.  3 clearly  shows  the  minimum  to  be  related  to  the  depletion  rear 
the  cloud  base  and  the  cloud  edge.  The  small  maximum  is  evidently  due 
mainly  to  the  convergence  of  vapor  evaporated  in  the  downdraft  and  the 
rain  shaft  and  carried  away  from  the  axis  by  the  shallow  outdraft  near 
the  ground  surface.  This  outdraft  is  shown  in  Fig.  13,  which  is  indi- 
cative of  the  circulations  in  Cloud  A even  though  the  data  are  actually 
from  Cloud  B.  There  is  also  a contribution  from  a little  "ledge"  that 
separates  the  subcloud  and  lateral  vapor  depletion  regions.  This  small 
anomaly  is  located  at  the  top  of  the  moist  layer  of  the  undisturbed 
atmosphere  and  in  a shearing  zone  between  the  mean  circulation  cells 
below  and  to  the  side  of  the  cloud. 

Brown  (1967),  following  the  suggestion  of  Braham  (1932)  defined 
the  available  energy  of  a cloud  as  the  latent  heat  that  has  been  re- 
leased by  condensation.  Using  a number  of  measurements  of  "cloud  mass" 
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(i.e.,  measured  cloud  volume  times  the  density  of  the  ambient  air)  and 
assumed  liquid  water  content  for  different  types  of  clouds,  he  developed 
an  empirical  formula  for  available  energy  as  a function  of  cloud  mass. 
Brown's  available  energy  is  comparable  to  the  change  of  latent  enthalpy 
computed  in  the  present  study.  However,  application  of  Brown's  formula 
to  Cloud  A yields  e value  of  14.5  TJ  as  compared  with  our  computed 
value  of  1.34  to  2.34  TJ  (see  Table  2).  Two  reasons  for  the  disparity 
may  be  mentioned.  Our  measure  of  "cloud  mass"  assumed  that  the  cloud 
was  a cylinder  of  600  ra  radius  and  4.6  km  depth,  which  is  the  approxi- 
mate envelope  of  the  cloud,  but  at  no  particular  time  during  the  com- 
putation did  the  cloud  entirely  fill  the  cylinder.  The  more  important 
reason  is  that  Brown  assumed  a uniform  and  rather  large  value  of  liquid 
water  content  throughout  the  cloud,  whereas  our  computations  yielded 
much  smaller  values  throughout  the  major  part  of  the  cloud.  In  view 
of  these  differences  and  of  the  fact  that  Brown  claimed  only  order-of- 
magnitude  accuracy  for  his  formula,  the  two  studies  working  from  en- 
tirely different  premises  give  comparable  results. 

THERMAL  ENTHALPY 

The  total  change  in  thermal  enthalpy  can  be  divided  into  two  terms, 
one  representing  water  vapor  (e.g.,  (3-18)),  the  other  dry  air  (e.g., 
(3-17)).  Since  no  liquid  was  present  at  t = 0 or  t = 60,  the  net  change 
in  thermal  enthalpy  of  liquid  was  zero.  The  change  in  the  density  of 
thermal  enthalpy  of  water  vapor  is  shown  in  Fig.  4.  The  pattern  largely 
reflects  the  change  in  water  vapor  content  rather  than  temperature 
changes,  and  therefore  is  similar  to  the  pattern  displayed  by  the  change 
in  latent  enthalpy.  The  size  of  the  change,  however,  is  about  one  order 
of  magnitude  smaller.  The  patterns  of  vertical  and  horizontal  movement 
of  thermal  enthalpy  of  water  vapor  are  shown  in  Figs.  8(b)  and  9(b). 
Again,  the  rearrangement  of  this  form  of  energy  is  similar  to  that  of 
latent  enthalpy. 

The  change  in  the  density  of  thermal  enthalpy  of  dry  air  is  shown 
in  Fig.  5.  Wien  is  used  in  energy  density  computations,  as  it  is  in 
all  categories  but  potential  energy,  the  mass  of  dry  air  remains  con- 
stant at  each  grid  point,  and  therefore  the  changes  in  enthalpy  of  dry 
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air  are  directly  related  to  changes  in  temperature.  It  is  noteworthy 
that  near  the  cloud  the  pattern  of  change  In  enthalpy  of  dry  air  is 
somewhat  similar  to  that  of  thermal  enthalpy  of  water  vapor  except  that 
the  sign  is  reversed.  The  cloud  evaporates  leaving  behind  wetter  but 
cooler  air.  The  magnitude  of  the  change  is  similar  (i.e.,  small  com- 
pared with  latent  enthalpy  changes). 

The  region  occupied  by  the  cloud  loses  energy  in  thp  form  of 
enthalpy  of  dry  air  (except  near  the  level  of  the  cloud  base)  and  the 
subcloud  region  (except  close  to  the  ground)  gains  energy  in  that  form. 
The  loss  of  energy  in  the  upper  part  of  the  cloud  directly  reflects  the 
cooling  by  evaporation  of  the  condensate  that  does  not  fall  out  as  rain. 
The  loss  of  energy  near  the  ground  is  caused  by  the  evaporative  cooling 
of  rain  in  the  downdraft  near  the  ground  and  the  outward  spread  of  this 
cooled  air.  A large  volume  of  air  receive < modest  increases  in  enthalpy. 
This  is  caused  by  the  (generally)  slowly  sinking  air  that  compensates 
for  the  upward  motions  within  the  cloud  itself.  The  region  of  warming 
includes  the  subcloud  layer,  except  near  the  ground,  and  extends  upward 
outside  the  cloud  to  a level  about  four-fifths  of  the  cloud  height  above 
the  cloud  base.  The  horizontal  extent  of  the  warming  is  at  least  5 
cloud  radii.  At  distances  greater  than  7 cloud  radii  (~  4200  m)  the 
pattern  becomes  less  distinct. 

The  vertical  rearrangement  of  enthalpy  of  dry  air  is  depicted  in 
Fig.  8(c),  which  clearly  shows  the  loss  of  energy  in  the  outflow  caused 
by  the  cool  downdraft,  the  gain  in  energy  in  the  air  between  the  down- 
draft  outflow  and  the  upper  mid -region  of  the  cloud,  and  tho  loss  of 
energy  in  the  vicinity  of  the  upper  region  of  the  cloud. 

The  horizontal  rearrangement  of  enthalpy  of  dry  air  is  shown  in 
Fig.  9(c).  Columns  of  air  from  the  axis  to  a point  somewhat  beyond 
the  cloud  radius  lose  energy,  and  an  increase  beyond  that  point  is 
clearly  shown.  The  fact  that  the  change  does  not  approach  zero  at  the 
outer  boundary  of  the  computational  domain  is  considered  an  unrealistic 
feature  related  to  the  closed  nature  of  the  domain.  The  closed  bound- 
aries presumably  lead  to  general  subsidence  that  is  greater  than  that 
found  in  nature.  That  is,  near  the  boundary  of  the  domain  the  air 
moves  downward  a greater  distance  than  would  occur  in  nature  and  the 
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effect  of  the  cloud  is  stronger  at  a greater  distance  than  would  be 
evident  in  nature. 

THERMAL  ENTHALPY  OF  MOIST  AIR 

It  was  noted  previously  t ha t changes  in  the  thermal  enthalpies  of 
dry  air  and  water  vapor  are  similar  but  of  opposite  sign.  The  water 
vapor  contribution  to  change  in  thermal  enthalpy  of  moist  air  will 
dominate  in  regions  close  to  the  cloud. 

The  vertical  rearrangement  of  thermal  enthalpy  of  moist  air  is 
shown  in  Fig.  8(d).  Tie  effect  of  evaporation  in  the  downdraft  below 
the  cloud  and  the  warming  outside  the  cloud  at  the  level  of  the  cloud 
base  are  clear  even  after  horizo  ital  summation.  The  increase  in  thermal 
enthalpy  in  the  slabs  at  the  height  of  the  tmper  portions  of  the  cloud 
is  caused  by  the  additional  water  content  injected  by  the  evaporation 
of  cloud  debris  at  that  level. 

The  horizontal  rearrangement  of  thermal  enthalpy  of  moist  air  is 
shown  in  Fig.  9(d).  The  contrasting  dominance  of  the  transport  of 
water  vapor  content  in  the  vicinity  of  the  cloud,  of  evaporative  cool- 
ing at  the  edge  of  the  cloud,  and  of  the  heatinp  caused  by  compensating 
sinking  motions  away  from  the  cloud  are  clearly  show. 

POTENTIAL  ENEF 

As  shown  in  rig.  6,  potertial  energy  increases  in  the  upper  portion 
of  the  cloud  and  out  to  about  5 cloud  radii,  and  decreases  be^  w the 
cloud.  ,.iere  are  decreases  in  potential  energy  over  most  or  the  domain 
external  to  the  cloud  except  directly  above  it.  The  absolute  value  of 
the  hange  in  energy  density  Is  not  a marked  function  of  radius,  except 
close  to  the  cloud  itself.  At  high  altitude,  even  near  the  outer  bound- 
ary of  the  domain,  absolute  values  of  potential  energy  change  are  of 

_3 

the  order  of  10  J m . This  is  a relatively  large  value  at  distances 
remote  from  the  cloud,  the  latent  enthalpy  and  thermal  enthalpy  changes 
being  about  one  order  of  magnitude  smaller. 

Potential  energy  density  is  the  only  category  computed  with  the 
use  of  an  air  density  that  was  allowed  ..o  vary  as  temperature  and  pres- 
sure vary  from  their  basic  states.  If  p rather  than  p were  used  in 
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computing  potential  energy  density,  as  it  is  for  other  energy  densities, 
there  would  be  little  variation  beyond  about  1.5  km  radial  distance,  as 
is  shown  by  Fig.  9(f).  However,  when  the  air  density  varies,  the  po- 
tential energy  varies  throughout  the  domain.  In  particular.  Fig.  9(e) 
shows  decreasing  potential  energy  with  increasing  radius,  illustrating 
the  influence  of  nearly  constant  subsidence  combined  with  the  effect  of 
the  square  of  radial  distance.  Figure  8(e),  using  p , and  Fig.  8(f), 
using  p^,  also  show  striking  differences  in  the  vertical.  In  general, 
with  variable  density,  as  is  required  for  consistency,  the  pattern  is 
smoother.  The  overall  change  in  potential  energy  is  negative,  the 
lower  portion  of  the  atmosphere  losing  more  energy  than  the  upper  por- 
tion gains. 

Some  further  questions  regarding  the  calculation  of  potential 
energy  are  discussed  in  Sec.  VI. 

STATIC  ENERGY 

The  overall  change  in  the  stability  of  the  atmosphere  is  indicated 
by  changes  in  the  static  energy,  which  is  the  sum  of  thermal  enthalpy, 
latent  enthalpy,  and  potential  energy,  'these  changes  are  shown  in 

7,  and,  not  surprisingly,  they  largely  reflect  changes  in  latent 
enthalpy.  The  lower  atmosphere  close  to  the  cloud  loses  energy,  while 
in  the  vicinity  of  the  upper  and  middle  portions  of  the  cloud  the  atmos- 
phere gains  energy.  This  situation  extends  horizontally  for  a distance 
of  several  cloud  radii  (and,  because  of  diffusion  and  mixing  with  the 
environment,  would  extend  even  farther  if  the  computation  continued  past 
the  lifetime  of  the  cloud).  This  rearrangement  of  energy  tends  to 
stabilize  the  atmosphere.  At  greater  distances  the  situation  is  re- 
versed. The  lower  atmosphere  receives  energy  while  the  upper  atmosphere 
loses  it.  This  causes  destabilization  of  the  atmosphere  and  increases 
the  likelihood  of  convective  activity.  The  magnitude  of  the  change 
toward  destabilization,  however,  is  nich  smaller  than  that  for  stabili- 
zation near  the  cloud.  Inspection  of  Fig.  / and  similar  figures  for 
changes  over  shorter  time  intervals  suggests  that  the  greatest  de- 
stabilization  occurs  at  a distance  of  about  5 cloud  radii  from  the 
axis  and  at  the  end  of  the  life  cycle  of  the  parent  cloud.  This  would 
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be  approximately  the  region  of  convergence  related  to  the  outward 
spreading  downdraft,  and,  consequently,  one  might  suppose  that  this 
would  be  a favorable  location  for  new  convective  activity.  The  dis- 
tance is,  of  course,  time-dependent  both  from  the  standpoint  of  changes 
in  vertical  energy  content  and  downdraft/outflow  convergence. 

Vertical  changes  in  static  energy  throughout  the  domain  (slab- 
summed)  are  shown  in  Fig.  8(g).  The  dominance  of  the  water  vapor 
transport  is  again  revealed  in  this  figure,  and  the  overall  stabili- 
zation of  the  atmosphere  is  apparent. 

Horizontal  changes  (tube-summed)  in  static  energy  are  shown  in 
Fig*  9(g).  Again,  the  dominance  of  water  vapor  transport  and  the 
latent  enthalpy  term  are  apparent.  The  cloud  is  clearly  shown  to  be 
a region  in  which  energy  is  concentrated  and  lifted  at  the  expense  of 
a large  surrounding  volume.  From  Fig.  9(g)  one  might  suspect  that  the 
effective  influence  of  the  cloud  extends  about  5 cloud  radii — a value 
similar  to  that  suggested  for  preferred  additional  convective  activity. 
Again,  this  effective  volume  is  time-dependent  and  the  present  remarks 
refer  to  the  situation  at  the  end  of  convection. 

ENVIRONMENT  STABILIZATION 

A key  function  of  cumulus  convection  is  to  relieve  potential  in- 
stability in  the  vertical  structure  of  the  atmosphere.  It  does  so  by 
transferring  static  energy  upward. 

The  static  energy  density  profile  (in  kilojoules  per  cubic  meter) 
of  the  environment  for  Cloud  A is  shown  in  Fig.  10,  together  with  the 
change  that  occurred  along  the  cloud  axis  as  a result  of  convective 
activity.  The  corresponding  changes  summed  over  entire  slabs  are 
shown  in  Fig.  8(g).  The  movement  of  static  energy  upward  is  apparent 
in  both  figures.  However,  if  the  specific  static  energy  (energy  per 
unit  mass  of  moist  air)  is  viewed,  as  in  Figs.  11  and  12,  a clearer 
picture  emerges. 

In  Fig.  11  a minimum  in  specitlc  static  energy  is  evident  at  a 
height  of  5 km.  Since  static  energy  is  a very  nearly  conservative 
property,  any  mass  of  air  rising  from  a level  lower  than  5 km  should 
remain  buoyant  until  it  reaches  at  least  5 km  regardless  of  mixing 
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Fig.  10-Static  energy  density  (lower  scale)  and  its  change  fnm  0 to  60  min 
(upper  scale)  on  the  central  axis  of  Cloud  A 

with  external  air  (entrainment).  If  no  entrainment  occurs,  the  parcel 
should  rise  until  it  reaches  a height  somewhere  above  5 km,  at  which 
the  specific  static  energy  equals  that  of  the  parcel.  If  entrainment 
occurs,  the  static  energy  of  the  rising  parcel  will  continuously  de- 
crease and  the  parcel  will  not  rise  so  high.  Strong  entrainment  may 
cause  the  static  energy  of  the  parcel  to  approach  that  of  the  environ- 
ment and  the  cloud  will  not  rise  much  above  the  level  of  minimum  static 
energy  (5  km  in  this  case).  If  the  cloud  is  poorly  organized  and  mix- 
ing and  diffusion  are  intense,  the  parcel  will  have  a specific  static 
energy  very  close  to  that  of  the  environment,  and  drag  forces  will 


Fig.  11— Static  energy  per  unit  mass  Fig.  12— Change  in  static  energy  per  unit  mass  of  moist  air  from  0 to  60  min  on  the  central  axis 

of  moist  air  at  time  0,  Clouo  A of  Cloud  A,  and  horizontally  averaged  over  the  entire  domain 
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p revert  it  from  rising.  Thus,  the  height  of  rise  of  the  cloud  tower 
is  dependent  on  the  degree  of  entrainment. 

Mixing,  or  entrainment,  is  a surface  effect,  whereas  the  total 
energy  in  a cloud  is  volume-dependent.  Thus,  the  larger  the  cloud 
radius,  the  less  will  be  the  mixing  in  p-oportion  to  total  energy  con- 
tent, and  the  higher  the  cloud  will  ’-ise.  In  some  cloud  models,  mainly 
those  in  one  dimension,  entrainment  is  explicit  (see,  e.g.,  Simpson 
°t  al-»  1965>-  In  these  models  the  height  of  rise  of  the  simulated 
cloud  is  controlled  by  specifying  the  cloud  radius,  and  consequently 
the  rate  of  entrainment.  In  the  cloud-simulation  model  used  here  en- 
trainment cannot  be  controlled  directly  because  it  is  implicit.  How- 
ever, the  radius  of  the  calculated  cloud  can  be  controlled  to  some 
extent  by  adjusting  the  horizontal  dimension  and  intensity  of  the 
perturbation  that  initiates  convection,  and  this  in  turn  implicitly 
affects  the  height  of  rise  of  the  cloud  top. 

In  the  present  experiment.  Cloud  A rose  to  j.S  km,  evidently  re- 
flecting moderate  entrainment.  Some  natural  clouds  in  the  area  at 
which  the  sounding  was  taken  (Puerto  Rico  on  20  August  1963)  were  of 
about  that  height,  but  a few  grew  larger,  and  many  more  smaller  ones 
were  also  present  (see  Figs.  5 and  7 of  Simpson  et  al.,  1965).  Thus, 
tht  simulation  must  be  considered  as  realistic,  at  least  as  regards 
cloud-top  height.  Whether  it  should  be  considered  as  a "median"  or 
"typical"  cloud  cell  with  regard  to  energetics  is  an  open  question, 
but  a reasonable  argument  could  be  made  to  that  effect. 

The  specific  static  energy  change  caused  by  the  development  of 
Cloud  A is  shown  in  Fig.  12.  Both  changes  along  the  cloud  axis  and 
throughout  the  domain  are  given.  The  upward  movement  of  energy  is 
clear,  particularly  in  the  region  of  active  cloud  growth.  The  large 
increases  in  the  upper  half  of  the  cloud  are  due  largely  to  increased 
latent  enthalpy.  At  these  levels  there  is  little  change  outside  the 
cloud  (Fig.  7),  but  there  is  an  average  increase  for  the  entire  domain 
because  of  the  in-cloud  contribution.  Near  the  level  of  the  cloud 
summit  but  away  from  the  cloud  there  is  some  increase  in  water  vapor, 
as  can  be  seen  in  Fig.  3;  but  outside  the  cloud,  at  the  height  of  its 
lower  part,  there  is  sufficient  loss  of  water  vapor  that  the  static 


I 


r 

-45- 

energy  change  for  the  domain  goes  negative,  although  it  is  still  posi- 
tive for  the  cloud  axis  (Fig.  12).  The  increase  in  vapor  near  the 
ground  several  cloud  radii  away  from  the  axis  causes  the  loss  of  static 
energy  in  the  entire  subcloud  layer  to  be  less  than  that  in  the  same 
layer  near  the  cloud  itself.  This  latter  loss  is,  again,  primarily 
due  to  loss  of  latent  enthalpy  (water  vapor).  Thus,  in  general  the 
latent  enthalpy  term  dominates  the  static  energy  change,  and  conse- 
quently it  dominates  the  con’ribution  of  the  convection  toward  stabiliz- 
ing the  atmosphere. 

There  are  other  important  effects,  however,  in  addition  to  changes 
in  latent  enthalpy.  In  particular,  the  potential  energy  increases  near 
the  cloud  and  decreases  in  the  outer  regions.  The  latter  changes  are 
particularly  important  in  relation  to  changes  in  jther  types  of  energy, 
which  are  much  smaller  in  that  region. 

It  is  clear  that  this  single  cloud  does  not  stabilize  its  environ- 
ment greatly,  although  the  tendency  is  certainly  in  that  direction. 
According  to  Fig.  11,  for  the  atmosphere  to  achieve  neutral  stability 

(static  energy  independent  of  height)  it  is  necessary  to  decrease  the 

4 -1 

static  energy  below  3 km  by  an  average  of  about  10  J kg  and  increase 

the  static  energy  above  that  level  by  a like  amount.  When  summed  over 

a cylinder  of  6 km  radius  and  3 km  depth,  with  appropriate  air  density, 

this  amounts  to  about  3 * 10^  J. 

In  the  present  experiment,  the  loss  sustained  between  the  ground 

surface  and  3 km,  determined  by  adding  slab-summed  differences,  was 
13 

only  1.2  * 10  J.  A cross-check  on  this  value  can  be  obtained  by 

noting  in  Table  2 that  the  total  condensate  of  3558  metric  tons  repre- 

6 6 12 

sents  a latent  enthalpy  of  3.558  * 10  * 2.5  * 10  , or  nearly  9 * 10  J. 

Since  these,  two  values  are  in  reasonable  agreement,  they  inspire  con- 
fidence in  the  validity  of  the  calculations.  It  follows  that  the 
energy  rearrangement  was  about  one-half  of  one  percent  of  that  required 
to  achieve  neutral  stability. 

If  this  experiment  had  shown  full  stabilization,  however,  it  would 
have  been  unrealistic,  for  in  nature  such  stabilization  is  achieved  not 
by  a single  cloud  tower  but  an  ensemble  of  clouds.  At  the  height  of 
its  vigor.  Cloud  A occupied  a horizontal  area  of  about  one  percent  of 
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ils  computnt ! onal  domain,  and  the  rest  of  the  area  was  cloud-free. 

This  is  a rattier  low  ratio  for  cloudy  area  to  clear  area,  as  many 
observations  show,  in  particular  those  of  Simpson  et  al.  (1965).  It 
was  mentioned  earlier  in  this  discussion  that  the  annulus  of  air  from 

3 km  to  6 km  radius  is  not  greatly  affected  by  the  cloud,  but  there 
is  a zone  ol  decreased  low-level  stability  at  a radius  of  about  i km. 

i b i s suggests  that  favorable  conditions  exist  lot  additional  convective 
activity  at  this  range.  Since  the  area  of  the  domain  beyond  3 km  is 

4 times  that  within  3 urn,  a case  can  be  made  that  the  area  controlled 
by  a convective  cell  of  this  size  is  about  one-fifth  the  area  of  this 
computational  domain,  so  that  in  nature  about  5 centers  of  convection 
would  occur  in  a domain  of  this  size.  The  geometry  of  the  model  pre- 
cludes the  appearance  of  these  aiditional  convective  centers  in  the 
simulation.  This  supposition  of  a natural  domain  does  not  conflict 
with  c_ nervations  of  cloud  density,  although  it  is  true  that  a tall, 
vigorous  cloud  can  depress  convection  in  its  immediate  vicinity. 

Another  consideration  is  that  Cloud  A is  not  really  representa- 
tive of  a natural  cloud  but  more  nearly  resembles  a cloud  tower.  Clouds 
of  this  depth  (5  km)  normally  consist  of  more  than  one  tower,  perhaps 
four  or  more.  Observations  from  aircraft  suggest  characteristic  active 
lifetimes  of  15  minutes  per  tower.  Thus  we  estimate  that  typically 
about  20  times  as  much  convection  would  be  occurring  in  this  domain  as 
is  simulated.  If  so,  then  in  one  hour  about  10  percent  of  the  total 
rearrangement  required  to  achieve  neutral  stability  would  be  accomp- 
lished. According  to  this  interpretation,  the  collective  effect  of 
convection  will,  because  of  the  interactions  between  clouds,  result  in 
a vertical  rearrangement  lying  somewhere  between  the  "on-axis"  and 
"domain"  patterns  shown  in  Fig.  12.  That  is,  the  entire  subcloud  layer 
will  lose  energy,  and  the  region  of  greatest  addition  of  energy  will 
be  somewhat  below  the  region  of  minimum  specific  energy. 

Clouds  with  a larger  volume-to-surf ace  ratio  than  the  model  cloud 
would  suffer  less  relative  mixing;  therefore  they  would  be  deeper  and 
more  energetic.  The  loss  of  energy  in  their  lower  regions  and  in  the 
subcloud  layer  would  also  be  greater.  Such  clouds  would  be  more  effec- 
tive than  the  cloud  modeled  in  transporting  energy  upward,  and  conversely 
smaller  clouds  would  be  less  effective. 
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The  deeper  clouds  are  required  if  energy  Is  to  be  transported  to 
the  higher  levels.  However,  under  normal  trade-wind  conditions  the 
larger  clouds  are  few  in  number,  and  it  is  not  obvious  whether  the  more 
numerous  but  individually  less  effective  small  clouds  or  the  rare  but 
individually  more  effective  large  clouds  contribute  more  to  the  ultimate 
rearrangement  of  energy  in  ’evels  of  the  atmosphere  occupied  by  both. 

It  should  be  noted  that  the  rearrangement  of  energy  and  the  degree 
of  stabilization  of  the  atmosphere  is  not  governed  by  the  amount  of 
rainfall  that  occurs.  A non-raining  cloud  will  still  cause  a rearrange- 
ment of  the  vertical  distribution  of  energy  in  the  atmosphere.  In  the 
model  cloud  the  amount  of  rainfall  is  an  order  of  magnitude  less  than 
the  total  condensation,  and  it  is  the  redistribution  of  water  vapor 
through  condensation,  advection,  and  evaporation  that  has  the  greatest 
effect.  In  the  case  of  Cloud  A,  the  enthalpy  lost  from  the  atmosphere 
by  the  fallout  of  rain  (thermal  enthalpy  of  water  removed  from  the 
system)  was  half  the  change  in  latent  enthalpy  and  much  less  than  the 
changes  of  either  thermal  enthalpy  of  the  atmosphere  or  potential 
energy;  see  Table  2.  The  product  of  the  mass  of  water  rained  out  and 
the  latent  heat  of  condensation  is  equal  to  the  loss  of  latent  enthalpy 
by  the  atmosphere. 
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V. ENERGY  TRANSPORT 


ENERCY_  FLUX 

The  flux  of  various  quantities  is  closely  related  to  the  wind 
field.  The  patterns  of  time-averaged  winds,  u and  w,  provide  a general 
picture  ot  the  horizontal  and  vertical  movement  of  the  air  and  conse- 
quently of  the  energy  transport.  These  fields,  for  Cloud  B,  are  shown 
in  Figs.  13  and  14.  The  strong  upward  movement  of  air  within  the 
cloud  (reaching  a maximum  of  nearly  1 m s_1)  and  the  strong  downward 
movement  beneath  the  cloud  (1.7  m s are  conspicuous.  The  mean  wind 
in  the  subcloud  downdraft  achieves  larger  values  than  the  mean  wind  in 
the  updraft  because  of  its  consistency  of  direction  rather  than  its 
maximum  instantaneous  strength. 

The  mean  horizontal  wind  pattern  shows  inflow  toward  the  cloud  in 
the  subcloud  layer  except  for  a region  of  strong  outflow  in  a shallow 
layer  near  the  ground  associated  with  the  deflection  of  the  downdraft 
by  the  rigid  boundary.  A maximum  of  nearly  1 m s'1  is  reached  in  this 
region.  Near  the  cloud  top  the  outflow  and  inflow  are  characteristic 
of  the  ring  vortex  which  the  cloud  circulation  to  some  extent  resembles. 
Air  movement  around  the  midsection  of  the  cloud  is  not  well  defined  in 
the  mean.  This  is  due  to  the  fact  that  air  at  a given  height  first 
moves  outward  as  the  cloud  top  approaches  from  below,  then  inward  as 
the  cloud  top  moves  through  the  zone;  consequently,  the  horizontal  wind 
components  tend  to  cancel  each  other  in  the  process  of  time  averaging, 
and  no  clear  pattern  of  mean  wind  flow  results. 

In  computing  energy  flux,  no  account  was  taken  of  the  changes  in 
the  form  of  energy  that  might  be  occurring  simultaneously.  Hence, 
conclusions  based  solely  on  flux  data  may  not  be  supported  by  data  dis- 
cussed in  Sec.  IV  regarding  the  overall  changes  in  enprgy. 

LATENT  ENTHALPY  FLUX 

The  flux  of  latent  enthalpy  averaged  over  a period  of  one  hour  is 
shown  in  Fig.  15.  The  strong  upward  movement  of  energy,  related  to 
moisture,  near  the  cloud  base  and  extending  through  the  cloud  volume, 
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is  clearly  shown.  In  the  cloud  region  the  size  of  the  term  w"nL"  (not 
illustrated)  is  about  one-third  the  total  average  flux,  whereas  away 
from  the  cloud  this  term  is  orders  of  magnitude  smaller  than  the  total 
average  flux.  Thus  the  correlation  between  the  departures  of  vertical 
wind  and  latent  enthalpy,  which  is  high  within  the  cloud  and  low  away 
from  the  cloud,  is  clearly  indicated.  The  horizontally  summed  average 
flux  of  latent  enthalpy  (AFL ) is  shown  in  Fig.  19(a).  There  is  vertical 
divergence  throughout  a layer  extending  from  near  the  ground  up  to 
1.8  km.  Above  the  height  of  the  cloud  base  (near  1.2  km)  latent  enthalpy 
flows  upward,  and  below  this  level  it  flows  downward. 

The  center  of  latent  enthalpy  moves  downward,  as  indicated  in 
Fig.  20,  suggesting  an  overall  downward  motion  of  water  vapor  through- 
out the  domain  over  the  lifetime  of  the  cloud.  The  general  subsidence 
in  the  outer  reaches  of  the  domain  and  the  rain  shaft  apparently  carry 
more  water  vapor  downward  than  the  active  convection  carries  upward. 

This  overall  movement  of  moisture  is  unexpected,  for  cumulus  convection 
normally  is  thought  to  cause  the  transport  of  moisture  upward.  This 
matter  was  further  examined  by  finding  the  height  at  which  half  the 
water  vapor  in  the  atmosphere  lay  above  and  half  below.  Initially, 
this  height  was  1310  m.  At  60  sec,  it  had  descended  a mean  distance 
of  10  cm  throughout  the  entire  domain,  but  on  the  cloud  axis  itself  it 
had  ascended  22  ra.  These  values  take  into  account  both  the  transforma- 
tion (to  rain)  and  the  advection  of  water  ' apor , and  they  corroborate 
the  flux  values.  There  is  reason  to  suspect  that  the  fixed  boundary 
and  the  geometry  of  the  simulation  is  causing  an  unrealistic  amount  of 
downward  flow  in  the  outer  region  of  the  domain,  thus  causing  the  net 
downward  flux  of  moisture.  Following  a similar  argument  developed  in 
the  discussion  of  environment  stabilization,  if  additional  clouds  were 
present,  the  overall  transport  probably  would  be  upward. 

THERMAL  ENTHALPY  FLUX 

The  flux  of  thermal  enthalpy  is  shown  in  Fig.  16.  The  pattern  and 
magnitude  of  the  Harare  similar  to  that  of  latent  enthalpy.  The  con- 
tribution of  the  w"r)g"  term  (not  illustrated)  is  about  1 percent  of  the 
total  flux  near  the  cloud  volume  and  diminishes  to  about  0.01 


percent 
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in  the  middle  of  the  eomput at ional  domain.  Thus,  unlike  the  latent 
enthalpy  flux,  there  is  little  correlation  between  the  departures  from 
the  mean  of  the  vertical  motion  and  thermal  enthalpy. 

Ihe  horizontally  summed  fluxes  of  thermal  enthalpy  of  dry  air, 
water  vapor,  and  condensed  liquid  are  shown,  respectively,  in  Figs. 
19(b),  19(c),  and  19(dl,  and  their  sum  is  shown  in  Fig.  19(e).  Flux 
for  the  gas  phases  changes  sign  at  the  level  of  the  cloud  base,  but  for 
the  liquid  phase  this  occurs  at  a somewhat  higher  level.  The  absolute 
value  of  the  flux  of  dry  air  enthalpy  is  an  order  of  magnitude  larger 
than  that  of  water  vapor  which,  in  turn,  is  an  order  of  magnitude 
larger  than  that  of  the  liquid  condensate. 

The  center  of  thermal  enthalpy  moves  downward  (Fig.  20). 

POTENTIAL  ENERGY  FLUX 

The  flux  of  potential  energy  is  shown  in  Fig.  17.  Again,  the 
pattern  and  magnitude  are  similar  to  those  of  latent  enthalpy,  but, 
as  in  the  case  of  thermal  enthalpy  flux,  there  is  relatively  little 
correlation  between  the  departures  of  vertical  wind  and  potential 
energy.  Immediately  aho-  e the  cloud  there  is  a slight  downward  flux 

in  potential  energy  occasioned  by  subsidence,  which  is  indicated  in 
Fig.  14. 

The  horizontal ly  averaged  potential  energy  flux  is  shown  in  Fig. 
19(f).  The  relatively  small  downward  flux  at  the  bottom  of  the  domain 
and  the  disproportionately  large  upward  flux  within  the  cloud  compared 
with  those  of  Intent  and  thermal  enthalpies  are  due  to  the  height  de- 
pendence of  potential  energy,  which  is  much  stronger  than  that  of 
enthalpy  and  generally  of  opposite  sign. 

The  center  of  potential  energy  moves  upward  (Fig.  20). 


STATIC  ENERGY  FLUX 

The  flux  of  static  energy  is  shown  in  Fig.  IS.  It  equals  the  sum 
of  the  potential  energy  and  thermal  and  latent  enthalpy  fluxes  and 
shares  their  characteristics.  Within  the  cloud  volume  is  about 

10  percent  of'wi^  whereas  elsewhere  in  the  domain  it  is  0.1  percent 
of  wms.  That  is  to  say,  vertical  wind  and  static  energy  are  slightly 
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correlated  in  the  region  of  maximum  convective  activity  and  hardly 
correlated  at  all  elsewhere. 

The  largest  component  of  the  static  energy  flux  is  the  thermal 
enthalpy  of  dry  air.  This  reflects  the  large  mass  of  air  involved  in 
the  flow  field  in  comparison  with  the  mass  of  water  substance.  The 
net  flux  of  dry  air  enthalpy  is  close  to  zero,  however,  and  the  up- 
ward movement  of  the  center  of  static  energy  (Fig.  20)  largely  reflects 
the  consistent  upward  flux  of  potential  energy. 

FLUX  OF  WATER  SUBSTANCE 

Since  some  of  the  major  forts  of  energy  depend  critically  on  the 
amount  of  water  vapor  or  liquid  coi densate  contained  in  the  air,  it  is 
instructive  to  see  how  these  constituents  are  transported.  Figure  21 
shows  the  vertical  flux  of  water  vapor,  and  Fig.  22  shows  that  of 
liquid  condensate.  The  familiar  pattern  of  the  mean  vertical  wind  is 
apparent  in  the  former.  The  latter  is  zero  over  most  of  the  domain 
because  of  the  absence  of  cloud  in  those  regions.  The  area  of  down- 
ward flux  above  5.5  km  and  to  the  side  of  the  cloud  is  associated  with 
very  small  amounts  of  water  carried  above  the  nominal  cloud  top  by  eddy 
diffusion,  which  enters  both  implicitly  through  the  differencing  algor- 
ithms and  explicitly. 

The  horizontally  summed  flux  values  for  water  vapor  and  total 
water  substance  (vapor  plus  condensate)  are  shown  in  Fig.  23.  The 
large  downward  flux  of  vapor  near  and  below  the  cloud  base  (1200  m)  is 
attributable  to  the  downdraft  beneath  the  cloud  as  wel]  as  subsidence 
in  the  environment  away  from  the  cloud.  The  predominance  of  the  vapor 
phase  with  reference  to  the  transport  of  water  substance  is  clearly 
evident  in  the  lower  part  of  the  domain,  but  near  the  cloud  top  the 
condensed  phase  accounts  for  about  half  of  the  flux.  Part  of  the 
reason  for  this  is  that  the  upward  flux  of  vapor  within  the  cloud  at 
a given  level  is  offset  to  a degree  by  the  downward  flux  outside  the 
cloud  at  the  same  level.  Also,  vapor  content  decreases  with  increas- 
ing altitude  throughout  the  domain,  but  over  much  of  the  region  occupied 
by  the  cloud  the  liquid  content  is  more  or  less  independent  of  height. 
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Fig.  23-Average  vertical  flux  summer!  by  slabs  (g  m'2  s'* ),  Cloud  B 

It  should  be  noted  that  the  flux  calculation  for  the  liquid  water 
does  not  take  into  account  the  gravitational  fall  of  the  larger  drops 
with  respect  to  their  associated  air  parcels. 

FLUX  DIVERGENCE 

The  vertical  divergence  of  several  categories  of  energy  may  be 
estimated  from  Fig.  19.  Wherever  a curve  in  that  figure  slopes  up- 
ward to  the  right  divergence  is  indicated,  and  wherever  it  slopes  up- 
ward to  the  left  convergence  is  indicated.  This  is  more  clearly  shown 
in  Fig.  24,  which  is  based  on  the  application  of  (3-41)  to  static 

energy.  In  effect,  the  curve  in  Fig.  24  is  the  derivative  of  the  curve 
in  Fig.  19(g). 


Vertical  energy  divergence  (kj  (200  m)"1  s-1) 

Fig.  24-Average  vertical  divergence  per  slab  of  static  energy  (KJ  (200  m)'1  s'1 ),  Cloud  B 


It  is  seen  that  static  energy  diverges  from  a region  extending 
from  near  the  ground  surface  to  a level  several  hundred  meters  above 
the  cloud  base.  The  level  of  maximum  divergence  is  at  cloud  base. 

At  higher  levels  divergence  alternates  with  convergence,  but  conver- 
gence predominates.  The  saw-tooth  pattern  of  alternating  positive 
and  negative  regions  of  divergence  in  Fig.  24  is  associated  with  a 
similar  variation  of  qv  with  height  in  the  original  sounding.  This 
is  apparent  in  the  dew-point  curve  of  Fig.  1.  This  irregularity  was 
the  reason  that  the  smoothed  version  of  the  sounding  was  used  for  much 
of  this  work. 


i 
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The  various  categories  of  enthalpy  and  the  potential  energy  that 
make  up  the  static  energy  all  have  patterns  of  divergence  similar  to 
that  shown  in  Fig.  24.  The  contribution  of  the  enthalpy  of  dry  air, 
however,  is  by  far  the  largest. 
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VI.  DIFFICULTIES— POTENTIAL  AND  REALIZED 


We  believe  that  some  comments  on  problems  that  often  occur  in 
studies  of  this  type,  or  that  occurred  during  the  course  of  this  work, 
would  be  valuable  in  guiding  future  efforts  and  provide  a feeling  for 
the  validity  of  these  results. 

ACCURACY  0F_  DIFFERENCES 

Energy  and  mass  differences  calculated  before  and  after  the  exist- 
ence of  a cloud  constitute  one  of  the  principal  quantitative  outputs. 
Often,  energy  contents  were  large  compared  with  differences.  This 
suggests  problems  concerning  the  significance  of  the  computed  differ- 
ences. A check  revealed  that  the  simple  subtraction  of  "before"  and 
"after"  values  of  all  quantities  was  accurate  to  at  least  two  decimal 
places.  This  is  probably  satisfactory  for  many  purposes,  but  not  for 
all;  consequently,  accuracy  was  substantially  improved  by  calculating 
differences  based  on  perturbations  from  a reference  state,  i.e., 

H = <J>T  - <fr0  - ($  + $'T)  - (7  + *'0)  - <fr'T  - $'0  - A*'  (6-1) 

Except  for  the  kinetic  energy  of  an  atmosphere  starting  at  rest,  the 
reference  value,  will  usually  be  large  compared  with  <j>'T  or  4>'q" 

Hence,  the  accuracy  of  the  differences  was  substantially  increased 
by  removing  the  large  reference  values  and  operating  with  the  residues. 
Another  device  used  to  increase  the  significance  of  the  differences 
was  the  factoring  out  of  large  constant  quantities  such  as  gz,  L,  or 
the  coefficients  of  specific  heat  (see  Sec.  III).  It  is  believed  that 
sufficient  care  was  thereby  taken  that  the  error  arising  from  taking 
small  differences  of  large  quantities  was  negligible. 

REARRANGEMENT  OF  ENERGY  CAUSED  BY  THE  FINITE-DIFFERENCE  SCHEME 

Normally,  the  cloud  model  has  been  run  using  real  atmospheric 
soundings;  consequently,  there  are  nonsmooth  transitions  in  the  vertical. 
The  finite-difference  scheme  used  in  the  computations  and  the  explicit 
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eddy  diffusion  terms  that  appear  in  the  model's  equations  both  tend 
to  smooth  out  these  sudden  transitions  of  water  vapor  content  or  temp- 
erature. Hence  there  is  a tendency  for  rearrangement  of  energy  (in 
this  case,  transport  of  thermal  enthalpy  and  latent  enthalpy  without 
transformation)  to  occur  regardless  of  the  presence  or  absence  of 
cumulus  convection.  This  smoothing,  which  can  be  significant  in  com- 
parison with  that  caused  by  convection,  appears  as  large  movements 
from  one  layer  to  an  adjacent  layer  across  the  entire  horizontal 
extent  of  the  domain.  It  is  distinguishable  from  convective  motion, 
which  will  have  a much  greater  variation  in  the  radial  direction  as 
well  as  a more  consistent  and  smooth  transition  in  the  vertical. 

Figure  25  illustrates  how  effects  caused  by  convection  can  be 
obscured  when  the  initial  sounding  is  not  smooth.  At  initial  time 
for  Cloud  B,  the  static  energy  summed  by  slabs  had  alternating 
shallow  layers  of  relatively  large  and  small  value  resulting  mainly 
from  large  vertical  gradients  of  humidity  in  the  original  unsmoothed 
sounding.  Subsequently,  numerical  smoothing  introduced  strong  nega- 
tive changes  at  levels  of  excess  relative  to  adjacent  grid  levels, 
and  strong  positive  changes  in  regions  of  deficit  relative  to  adjacent 
grid  levels,  resulting  in  the  saw-toothed  curves  of  Fig.  25.  This 
problem  may  be  reduced  by  using  a smoothed  sounding,  as  was  done  for 
Cloud  A. 

Compai ison  of  Fig.  8(g)  for  Cloud  A with  Fig.  25  for  Cloud  B shows 
that  the  use  of  a smoothed  sounding  not  only  produces  a smoother  curve 
of  slab-summed  energy  change  versus  height,  but  also  reduces  the  change 
for  a given  slab  by  an  order  of  magnitude  or  more.  In  this  connection, 
however , it  should  be  noted  that  even  the  large  changes  of  Fig.  25 
represent  only  1 percent  of  the  basic  values.  Furthermore,  the  posi- 
tive and  negative  changes  tend  to  counterbalance  each  other,  so  the 
total  change  cannot  be  too  far  off  the  mark,  at  least  in  a relative 
sense.  Nevertheless,  the  results  with  the  smoothed  sounding  inspire 
greater  confidence. 

REARRANGEMENT  OF  ENERGY  CAUSED  BY  EDDY  DIFFUS ION 

Usually,  cloud  dynamics  models  employ  an  eddy  diffusion  term  in 
most,  or  all  of  the  prognostic  equations.  Its  purpose  is  to  account  for 
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Fig  25— Change  from  0 to  15  min,  0 to  30  mm,  0 to  45  min,  and  0 to  60  min  in  static  energy 
summed  by  slabs,  from  the  unsmoothed  sounding  used  in  Cloud  B 


son  les  of  motion  too  small  to  be  resolved  by  the  finite  grid.  In  this 
model  the  eddy  diffusion  terms  have  the  form  where  V is  the  eddy 

diffusion  coefficient  and  •{’  might  be  the  momentum,  temperature  depar- 
t ire,  or  mixing  ratio  of  a constituent.  The  effect  of  eddy  diffusion, 
at  least  in  this  form,  is  to  "smear"  values  even  in  the  absence  of  a 
wind  field. 
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Whether  or  not  eddy  diffusion  is  to  be  used  is  a matter  of  judgment 
that  ultimately  depends  on  whether  the  cloud  simulation  is  more  real- 
istic with  or  without  the  term.  We  have  not  finally  resolved  this 
matter,  but  our  best  results  seem  to  have  come  with  the  use  oi  a coef- 
ficient of  v MO  id  sec  1 applied  to  the  momentum,  temperature  depar- 
ture, mixing  ratio  of  water  vapor,  and  mixing  ratio  of  the  part  of 
liquid  that  is  in  the  form  of  small  droplets  and  V = 0 to  the  part  of 
the  liquid  that  is  in  the  form  of  large  drops. 

The  model  starts  its  computations  with  a specified  atmosphere  that 
is  uniform  in  the  horizontal.  To  induce  convection  some  kind  of  per- 
turbation must  be  introduced.  The  perturbation  might  be  a horizontal 
gradient  in  the  buoyancy  in  some  limited  region,  which  will  induce 
motion,  or  it  might  be  a field  of  motion  itself.  In  both  cases  self- 
sustaining  convection  can  be  set  off  only  by  updraft  sufficient  to 
cause  some  condensation. 

Normally,  a buoyancy  perturbation  arising  from  altered  humidity 
has  been  employed  in  using  this  model.  It  is  a reasonable  emulation 
of  nature  s methods  of  inducing  convection  (see  Bunker  et  al.,  1949, 
and  Malkus , 1954),  and  it  was  used  in  the  computation  of  Cloud  B.  It 
is  equally  reasonable  to  start  with  an  initial  field  of  motion;  this 
was  done  in  Cloud  A. 

Any  perturbation  amounts  to  a perturbation  in  energy.  Because 
the  initial  state  of  the  atmosphere  includes  the  perturbation  insofar 
as  the  energy  calculations  are  concerned,  some  of  the  energy  changes 
can  be  affected  by  the  spread  of  that  perturbation  through  numerical 
diffusion  from  the  finite-difference  scheme  or  through  eddy  diffusion. 
Even  though  the  eddy  diffusion  terms  represent  real  phenomena,  their 
finite-difference  implementation  can  cause  the  spreading  of  small 
quantities  of  the  substance  involved  at  speeds  of  one  grid  unit  per 
time  step,  well  in  excess  of  any  speed  actually  permitted  by  the 
model.  This  spreading,  which  is  more  apparent  with  a humidity  per- 
turbation than  with  a field  of  motion,  can  occur  even  in  the  absence 
of  convection.  And,  as  in  the  case  of  vertical  irregularities  in  the 
basic  sounding,  the  energy  rearrangements  occasioned  by  the  spreading 
oi  the  perturbation  may  be  comparable  in  size  to  those  caused  by  the 
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convection.  Aga in , the  problem  is  less  serious  with  respect  to  flux 
calculat ions. 

Ihe  phenomenon  is  illustrated  in  Fig.  26.  The  initial  humidity 
perturbation  had  a maximum  at  the  central  axis  and  decreased  to  zero 
at  2 km  radial  distance.  Artificial  diffusion  causes  an  unusually 
large  transport  of  humidity  outward  across  the  edge  of  the  perturba- 
tion, resulting  in  a cusp  at  2 km  in  all  the  curves  of  Fig.  26. 

To  minimize  this  effect  the  energy  of  the  perturbation  should  be 
minimal.  The  humidity  perturbation  contains  large  amounts  of  energy 
in  the  form  of  thermal  and  latent  enthalpy;  hence  it  is  not  optimum 
for  this  purpose.  In  the  present  experiment  at  least,  kinetic  energy 
is  one  of  the  smallest  forms  of  energy  in  the  atmosphere,  and  further- 
more, if  free  convection  is  not  initiated,  kinetic  energy  is  rapidly 
damped.  It  does  affect  the  flux,  however. 

Lt  was  as  a consequence  of  these  considerations  that  Cloud  A with 
a momentum  impulse  was  used  to  analyze  the  energy  changes  and  Cloud  B 
with  humidity  impulse  was  used  for  the  flux  calculations . 

COMPUTATION  GEOMETRY 

The  computational  domain  is  a cylinder  with  a vertical  axis.  It 
is  described  by  cylindrical  coordinates  with  a uniform  mesh  length  D 
in  the  horizontal  and  the  vertical.  Axial  symmetry  is  assumed,  so 
each  interior  grid  point  represents  an  annulus  of  width  and  depth  D 
and  radius  depending  on  the  horizontal  index  of  the  grid  point.  (Eound- 
ary  points  are  special  cases.)  It  follows  that  the  volume  represented 
by  a grid  point  varies  with  the  square  of  the  radial  distance,  and  that 
small  differences  in  energy  density  at  large  radial  distances  can 
assume  great  importance  when  summing. 

Ideally,  the  computational  domain  should  be  so  large  that  the 
atmosphere  at  the  outer  and  upper  boundaries  is  undisturbed  bv  the  con- 
vection. In  this  work  the  outer  boundary  is  at  6 to  8 cloud  radii, 
assuring  relatively  little  disturbance,  but  permitting  some.  The  in- 
creases noted  toward  the  outer  boundary  in  some  of  the  graohs  of  Fig.  9 
are  probably  the  result  of  rather  insignificant  changes  magnified  by 
the  effect  of  the  square  of  the  radial  distance.  Caution  must  be  used 
in  assessing  the  realism  and  significance  of  some  of  the  results. 
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Fig  26  Change  from  0 to  10  min,  0 to  20  min,  0 to  30  min,  0 to  40  min,  0 to  50  min,  and  0 to 
60  min  of  static  energy  summed  by  tubes,  from  a run  with  a humidity  perturbation 
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The  geoir.et.ry  and  resolution  are  such  that  the  cloud  top  and  side 
are  well  defined  by  the  energy  differences.  All  values  show  a sharp 
decrease  across  the  cloud  bounda.  es. 

POTENTIAL  ENERGY— AIR  DENSITY  IN  A BOUSSINESQ  MODEL 

Since  the  cloud  simulations  on  which  this  study  of  energetics  is 
based  come  from  a Eoussinesq  model,  the  expressions  for  the  various 
types  of  energy  should  be  formulated  in  a way  that  is  consistent  with 
the  Boussinesq  assumptions.  In  the  present  study  this  was  done  by 
deriving  the  energy  equations  (2-32)  for  unit  mass  of  cloudy  air  and 
(2-40)  for  unit  volume  from  the  Boussinesq  eeuation  of  motion.  The 
outcome  was  that  wherever  air  density  enters  iu  appears  as  a reference 
density  dependent  only  on  height,  except  that  in  the  potential  energy 
term  the  full  current  local  density  appears.  In  all  cases  the  air 
density  is  dependent  on  the  water  vapor  and  condensed  water  content, 
but  in  the  case  of  potential  energy  it  is  also  dependent  on  pressure 
and  temperature. 

To  evaluate  (2-42),  one  must  know  the  value  of  p,  which  in  accord- 
ance with  (2-10)  is  equal  to  p + p' . The  reference  density  is  computed 
as  a function  of  z once  and  for  all  at  initial  time.  The  density  de- 
parture can  be  found  with  the  use  of  the  Boussinesq  equation  of  state 
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The  cloud  model  does  not  in  principle  require  that  pT  be  computed, 
for  as  a rule  p'/p  is  at  least  an  order  of  magnitude  smaller  than  T'/T. 
For  some  purposes,  however,  such  as  the  present  one,  it  is  desirable 
to  know  the  value  of  p'.  This  can  be  found  by  solving  an  elliptical 
partial  differential  diagnostic  equation,  or  more  simply,  by  solving 
the  third  equation  of  motion  for  Dp'/rlz  with  dw/dt  a known  quantity 
from  values  of  w at  the  beginning  and  end  of  the  time-step,  followed 
by  integrating  from  the  rigid  lid  of  the  computational  domain  (at 
which  p'  = 0)  to  the  ground  surface.  The  disadvantage  of  solving  for 
pressure  in  this  way  is  that  the  solution  lags  one  time-step  behind 
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the  microphysical  calculations,  and  so  there  is  a potential  for  con- 
siderable error  Since  p'  is  so  small  anyhow,  the  error  should  not 
be  troublesome. 

In  order  to  determine  the  importance  of  the  complete  air  density 
value  to  the  potential  energy  density,  some  calculations  were  made 
using  pa  rather  than  in  (2-42).  This  is  called  a "constant  density" 
situation  because  air  density  does  not  depend  on  the  current  tempera- 
ture and  pressure  even  though  it  varies  with  height  and  vapor  and 
liquid  content.  As  can  be  seen  from  Fig.  27,  the  pattern  of  change  of 
potential  energy  using  "constant"  density  closely  resembles  that  of 
latent  enthalpy  (Fig.  3),  bit  the  values  are  several  orders  of  magnitude 
smaller.  This  would  lead  one  to  conclude  that  changes  in  potential 
energy  are  small  and  relatively  unimportant. 

The  results  are  significantly  different,  however,  if  potential 
energy  density  is  computed  with  p in  (2-42),  as  it  should  be,  and 
with  the  use  of  (6-2);  see  Fig.  6.  The  absolute  value  of  potential 
energy  change  is  an  order  of  magnitude  larger  than  in  Fig.  27,  and  al- 
though each  figure  shows  a maximum  in  the  vicinity  of  the  cloud,  they 
differ  elsewhere. 

F.ven  though  the  use  of  variable  dry-air  density  in  computing  po- 
tential energy  is  a more  consistent  application  of  the  Boussinesq 
approximation  than  is  the  use  of  constant  dry-air  density,  it  still  is 
not  entirely  satisfactory  in  the  present  formulation.  Although  the 
atmosphere  should  suffer  a net  loss  in  energy  due  to  the  removal  of 
enthalpy  of  rain  water,  the  loss  shown  by  the  present  computation  is 
seven  times  too  large.  This  loss  appears  to  be  related  to  large  de- 
creases in  potential  energy  occurring  near  the  outer  edge  of  the  com- 
putational domain.  This,  in  turn,  is  probably  related  to  widespread 
subsidence  resulting  from  continuity  in  a closed  domain  of  arbitrary 
size.  It  is  conceivable  that  the  problem  might  be  lessened  if  the 
model  permitted  the  spacing  between  neighboring  convective  cells  to 
be  determined  dynamically. 

Also,  when  using  variable  dry— air  density  there  is  no  assurance 
that  the  mass  of  dry  air  is  conserved.  If  there  is  a loss  in  mass, 
there  will  be  a loss  in  potential  energy.  This  could  account  for  the 


NOTE  Shaded  areas  represent  a decrease. 


Fig  27— Change  from  0 to  60  min  of  potential  energy  density  (J  m ' hr  1 ) 
using  p a(qa  + qv  + <i^)  for  total  air  density,  Cloud  A 


apparently  too  large  decrease  in  the  potential  energy  of  the  domain. 
The  large,  region  of  Fig.  6 showing  decreased  potential  energy  density 
with  little  radius  dependence  suggests  an  overall  loss  of  mass.  Per- 
haps dry-air  mass  conservation  should  be  accounted  for  in  a manner 
similar  to  that  used  to  conserve  water  substance;  t tis  probably  would 
improve  the  energy  balarce. 
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VII.  CONCLUS IONS 


Although  results  from  only  two  cloud  runs-  are  shown  in  the  analysis, 
data  derived  from  runs  using  other  soundings  or  variations  in  some  de- 
tails of  the  model  were  available  to  supplement  these  data.  Also,  the 
numerical  results  generally  support  (i.e.,  compare  favorably  with)  the 
analysis  of  Braham  (1952)  based  on  field  observations.  We  therefore 

feel  that  °ur  ^sults  warrant  general  interpretation— with  reservations, 
of  courne. 

This  study  involves  the  energetics  of  an  isolated  cloud  tower. 

Most  convection  occurs  with  some  degree  of  mesoscale  (or  even  synoptic- 
scale)  organization,  and  cells  do  not  occur  in  isolation.  How  the 
transformation  and  rearrangement  of  energy  is  affected  by  the  inter- 
ference of  superimposed  effects  is  difficult  to  judge.  The  nature  of 
the  computational  domain,  a closed  box,  also  limits  the  generalization 
that  may  be  drawn.  The  effect  of  the  domain  is  most  apparent  near  the 
outer  boundary,  and  unfortunately  it  is  here  that  large  masses  of  air 
are  represented. 

Another,  perhaps  more  important,  source  of  uncertainty  in  apply- 
ing these  results  is  the  lack  of  balance  in  the  energy  summation.  While 
the  computations  of  total  energy  before  and  after  convection  agree  to 
one  part  in  10  , the  fact  remains  that  significant  changes  in  the  re- 
sults could  come  about  from  changes  in  the  distribution  and  form  of 
that  one  part.  We  believe  that  if  the  computation  were  exactly  con- 
servative, the  patterns  of  change  would  essentially  duplicate  those  of 
this  study.  However,  the  vertical  and  horizontal  summations  might  show 
important  differences.  Evidence  suggests  that  the  most  questionable 
terms  of  the  present  computation  are  the  potential  energy  and  the  ther- 
mal enthalpy  of  dry  air. 

FINDINGS  PERTINENT  TO  CUMULUS  PARAMETERIZATION 
IN  MODELS  OF  LARGER  SCALE  ~ “ 

A finding  of  particular  importance  to  the  problem  of  parameteriza- 
tion of  cumulus  convection  is  that  the  energy  available  to  stabilize 
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the  atmosphere  is  equal  to  the  total  latent  heat  released  during  con- 
densation (to  generalize,  include  deposition)  and  not  to  the  latent 
energy  lost  through  removal  of  water  from  the  atmosphere  in  the  con- 
densed phase  as  rain  or  snow.  This  confirms  Braham's  (1952)  conclu- 
sions. He  showed  from  data  of  the  Thunderstorm  Project  that  the  latent 
heat  of  condensation  of  the  liquid  remnants  of  "he  cloud  plus  that  of 
condensate  that  was  re-evaporated  considerably  exceeds  that  of  the  rain 
that  fell  to  the  ground.  Stated  another  way,  it  is  hazardous  to  esti- 
mate (or  parameterize)  the  amount  of  rainfall  from  the  energy  required 
to  stabilize  the  atmosphere.  Such  an  approach  is  basic  to  the  general 
circulation  model  that  has  been  used  by  Rand  to  model  climate  (Gates 
et  al.,  1971).  Experience  (Gates,  1972;  Rapp  and  Warshaw,  1974)  has 
shown  that  this  model  overestimates  the  amount  of  rainfall.  From  our 
analysis  of  cloud  energetics,  such  an  overestimate  is  to  be  expected. 

In  effect,  the  rainfall  parameterization  assumes  the  cloud  to  have  a 
precipitation  efficiency  of  100  percent  (rain  production  equals  total 
amount  of  condensation).  By  contrast,  in  the  present  calculation  the 
precipitation  efficiency  was  only  8 percent. 

The  present  work,  in  agreement  with  Brahan,  suggests  that  the 
equivalent  amount  of  condensate  might  be  taken  as  the  amount  of  heat 
required  to  stabilize  the  atmosphere  and  that  the  amount  of  rainfall 
should  be  a fraction  of  the  total  condensate  (i.e.,  the  condensate 
multiplied  by  the  "cloud  efficiency").  Cloud  efficiency  will  vary 
with  circumstances.  The  larger  the  cloud,  in  general,  the  greater  the 
efficiency.  Nominal  values  run  between  0.25  and  0.75  approximately, 
but  in  certain  climatic  regions  (e.g. , the  Southwestern  United  States) 
even  large  storms  are  observed  not  to  produce  rain  at  the  ground  be- 
cause of  substantial  evaporation  in  the  deep,  dry  subcloud  layer. 

Results  of  a run  that  used  the  initial  conditions  of  Cloud  A but 
did  not  permit  precipitating  drops  to  form  suggest  that  nonraining 
clouds  do  not  stabilize  the  subcloud  region  of  the  atmosphere.  In- 
stead, subsidence  and  downward  transport  of  water  vapor  tend  to  in- 
crease the  static  energy  below  the  cloud.  The  cloud  layer  itself, 
however,  is  stabilized.  If  this  is  true  in  nature,  one  would  expect 
that  small  convective  elements,  probably  "bubbles"  of  the  type  en- 
visioned by  Scorer  and  I.udlam  (1953),  would  effect  the  stabilization 


of  the  subcloud  region.  These  convective  elements  would  not  necessarily 
result  In  condensation  and  would  probably  be  below  the  limits  of  resolu- 
tion of  the  present  model.  A cumulus  model  incorporating  a good  boundary- 
layer  parameterization  might  demonstrate  this  process. 

The  present  simulations  suggest  that  the  rate  at  which  convective 
processes  effect  the  adjustment  of  the  atmosphere  toward  stability  is 
slow,  a period  of  the  order  of  ten  hours  being  required  to  decrease  the 
instability  to  1/e  of  its  original  value.  This  is  in  sharp  contrast  to 
the  value  of  one  hour  used  in  the  original  version  of  the  Rand  climate 
model  (Gates  et  al. , 1971).  The  indication  by  the  present  cloud-model 
calculations  that  the  rate  of  adjustment  is  slow  is  in  accord  both  with 
general  observations  of  convective  activity  and  with  Braham's  (1952) 
conclusions. 

FINDINGS  REGARDING  THE  RELATIVE  IMPORTANCE  OF  THE  ENERGY  TERMS 

The  kinetic  energy  involved  in  cumulus  convection  is  negligible  in 
comparison  with  static  energy.  The  largest  component  in  the  static 
energy  is  the  enthalpy  of  dry  air.  The  latent  enthalpy  is  about  one 
order  of  magnitude  smaller,  but,  because  of  the  stratification  of  water 
vapor  in  the  atmosphere  and  the  involvement  of  water  vapor  in  the  cloud- 
forming  process,  the  latent  enthalpy  term  is  much  more  variable  in  time 
and  space  and  therefore  accounts  for  the  majority  of  the  observed  changes 
in  energy.  The  water-vapor  thermal  enthalpy  term  generally  varies  with 
latent  enthalpy,  but  is  about  one-fifth  its  value. 

The  potential  energy  term  is  important,  particularly  in  locations 
away  from  the  immediate  vicinity  of  the  cloud.  The  overall  loss  of 
potential  energy  is  related  to  the  large  volume  of  subsiding  air  in  the 
outer  part  of  the  domain.  If  there  had  been  five  convective  elements 
in  a domain  of  this  size  (a  reasonable  number,  as  has  been  previously 
indicated),  this  loss  presumably  would  have  been  much  diminished. 

FINDINGS  REGARDING  THE  LIFE  CYCLE  OF  CONVECTIVE  ELEMENTS 

During  the  active,  growing  period  of  the  cloud  essentially  all  the 
condensate  that  is  to  form  is  produced.  Thus  the  energy  used  for  atmos- 
pheric stabilization  becomes  available  during  this  time.  Since  the  end 
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of  the  growing  phase  signals  the  end  of  new  energy  production,  a rapid 
decline  in  the  energy  available  for  rearrangement  is  experienced  as 
evaporation  continues.  In  the  runs  used  for  this  study,  the  cloud 
reaches  maturity,  and  nearly  all  of  the  significant  changes  in  energy 
forms  occur  during  the  first  30  minutes  of  simulated  time.  During  the 
subsequent  period  of  decay  there  is  some  rearrangement  of  energy  due 
to  the  field  of  motion,  evaporation,  and  numerical  smoothing,  but  it 
is  of  minor  importance. 
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